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ABSTRACT 


ital  Multivariable  Control  Laws  For  The  KC-135A 


I 

,  \ 

1  Multivariable  design  technique^ developed  by  Professor 
Brian  Porter  of  the  University^of  Salford,  England,  are  used 
to  develop  digital  control  laws  for  the  KC-135A.  Control 
laws  are  developed  for  each  of  three  diverse  flight 
conditions.  MULTI,  a  computer  program  used  in  the  design, 
is  modified  to  account  for  computational  tim^  delay.  The 
effects  of  a  computational  time  delay  on  the  controllers 
developed  are  presented. 

The  controllers  developed,  using  these  techniques,  ut-^-C 

> 

‘utilize1  output  feedback  with  proportional  plus  integral 
control.  Because  of  the  structure  of  the  system, 
measurement  variables,  in  addition  to  the  outputs,  are 
necessary.  A  reduced  controller  is  modeled  by  setting  some 
of  the  required  feedback  gains  to  zero.  A  comparison  of  the 
results  of  the  reduced  controller  to  the  complete  controller 
is  presented. 

A  robustness  controller  is  tested  by  performing  specific 
maneuvers  at  more  than  one  flight  condition.  The  robust 
controller  is  of  the  reduced  form  and  is  required  to  perform 
under  the  constraints  of  the  computational  time  delay. 

This  thesis  concludes  that  actuator  dynamics  play  a 

i 

significant  role  in  the  development  of  control  laws,  and, as 


a  result, 2nd  order  or  higher  actuator  models  should  be  used 


in  future  studies. 


DIGITAL  MULTIVARIABLE  TRACKER  CONTROL 
LAWS  FOR  THE  KC-135A 

I .  INTRODUCTION 


Background 

Aircraft  are  now  being  designed  with  closed- loop 
fly-by-wire  flight  control  systems.  These  multiple 
input/multiple  output  (MIMO)  systems  are  required  to  achieve 
stability  and  desirable  performance  while  transforming  the 
pilot's  command  inputs  into  movement  of  the  flight  control 
surfaces.  At  the  heart  of  this  MIMO  system  is  a  digital 
computer . 

Presently,  controllers  for  such  MIMO  systems  are  often 
designed  using  classical  single  input/single  output  (SISO) 
and  Linear  Quadratic  Gaussian  (LQG)  methods.  Each  of  these 
methods  have  some  serious  drawbacks.  The  well  known 
classical  SISO  method  may  yield  a  satisfactory  design, 
however,  it  allows  for  examination  of  only  a  single 
input/output  combination  at  any  one  time.  LQG  methods  can 
handle  the  MIMO  case,  but  all  internal  states  have  to  be 
accessible  for  the  design  method  to  work.  If  not 
accessible,  LQG  requires  the  use  of  complicated  state 
estimators  to  provide  the  system  with  an  estimate  of  the 
unaccessible  states. 

Professor  Brian  Porter  and  his  associates  of  the 
University  of  Salford,  England,  have  suggested  the  use  of 
four  direct  design  methods  capable  of  synthesizing  a 
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controller  which  can  handle  MIMO  systems  while  at  the 

* 

same  time  providing  for  a  high  degree  of  output  decoupling  ■< 

(Ref.  2).  These  design  techniques  produce  controllers  for  5 

highly-interactive  linear  multivariable  systems.  ^ 

) 

Initial  attempts  to  synthesize  an  aircraft  flight  a 

3 

control  system  using  Porter's  techniques  has  proven  to  be 
successful  (Ref.  3  and  4).  These  investigations  found  the 
techniques  to  be  fairly  straightforward.  Still,  much 
remains  to  be  studied  before  a  final  design  can  be 
implemented  on  future  aircraft. 

AFIT  students  have  developed  a  computer  program  called 
MULTI  (Ref.  4)  that  utilizes  the  design  techniques  of  Brian 
Porter.  MULTI  allows  the  user  to  both  design  and  simulate  a 
MIMO  controller  with  the  help  of  a  digital  computer.  The 
computer  provides  a  means  for  multiple  iterations  of  the 
design  and  simulation  to  allow  the  user  to  achieve  a 
suitable  design. 

Problem 

The  purpose  of  this  thesis  is  to  investigate  the 
development  of  tracker  control  laws  for  the  KC-135A 
aircraft.  Controller  performance  is  demonstrated  through 
the  use  of  the  computer  program  MULTI  by  performing 
designated  flight  maneuvers.  Controllers  are  developed  for 
each  of  three  flight  conditions.  For  each  flight  condition, 
four  maneuvers,  two  lateral  and  two  longitudinal  are 


o 


performed.  A  robust  controller  is  developed  which  is 
capable  of  producing  satisfactory  responses  for  a  single 
maneuver  for  all  three  flight  conditions  studied. 

This  thesis  also  incorporates  a  modification  to  MULTI 
that  allows  for  the  addition  of  a  computational  time  delay 
when  simulating  the  results  of  a  controller.  Controller 
performance  is  tested  against  this  computational  time  delay. 

Approach 
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In  order  to  accomplish  the  desired  results,  six  steps 
are  taken  in  completion  of  this  thesis. 

1.  Development  of  the  aircraft  equations  of  motion  for 
three  separate  flight  conditions,  including  a  body 
bending  mode  in  the  longitudinal  case. 

2.  Adding  to  MULTI  the  capability  to  include,  at  the 
user's  request,  a  computational  time  delay  when 
simulating  design  results. 

3.  Design  of  digital  multivariable  tracker  control 
laws  using  the  appropriate  design  methods  of 
Chapter  II. 

4.  Evaluation  of  the  effects  of  computational  time 
delay  on  each  of  the  controllers  developed. 

5.  Investigation  of  minimum  gain  designs  where  the 
gain  Matrices  Kg  and  K-^  are  minimized  for 

each  maneuver.  L 

6.  Development  of  a  robust  controller  capable  of 
performing  a  single  maneuver  for  all  flight 
conditions. 


Assumpt ions 

It  is  assumed  that  computer  simulations  give  ; 
realistic  approximation  of  aircraft  motion.  Other 
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assumptions  as  necessary  in  the  development  of  certain  parts 
of  this  thesis  are  listed  as  they  are  needed. 

Presentat  ion 

Six  chapters  are  included  in  this  thesis.  Chapter  II 
describes  the  theory  underlying  the  design  method  used  in 
this  study.  Chapter  III  describes  the  KC-135A  aircraft  and 
Appendix  A  details  the  development  of  the  necessary  aircraft 
equations  of  motion.  Chapter  TV  present  the  computational 
time  delay  modification  made  to  MULTI.  Chapter  V  describes 
the  development  of  tracker  control  laws  for  the  KC-135A  for 
the  three  different  flight  conditions.  This  includes 
discussions  of  the  results  obtained  for  flight  condition  #3 
using  computational  time  delays  and  minimum  K  designs.  Also 
included  in  Chapter  V  is  a  presentation  of  the  robust  con¬ 
troller  developed  for  each  of  the  four  maneuvers.  Chapter 
VI  is  the  conclusion  and  offers  recommendations  for  future 
studies  and  suggestions  for  possible  improvements  and/or 
additions  to  MULTI. 
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II.  MULTIVARIABLE  CONTROL  LAW  THEORY 

Introduction 

The  control  law  theory  used  in  this  thesis  was 
developed  primarily  by  Professor  Brian  Porter,  of  the 
University  of  Salford,  England  with  the  help  of  Doctor  A. 
Bradshaw,  also  of  the  University  of  Salford.  Professor 
Porter's  research  was  performed  under  a  contract  sponsored 
by  the  Air  Force  Flight  Dynamics  Laboratory  (AFWAL/FIGL) , 

Wr ight-Patter son  Air  Force  Base,  Ohio.  Only  a  summary  ot 
the  theory  on  singular  perturbation  multivariable  controller 
design  is  presented  in  this  chapter  and  in  Reference  2.  A 
truly  interested  reader  can  trace  the  development  of  this 
theory  from  its  early  stages,  by  researching  the 
bibliography  in  chronological  order. 

The  multivariable  controller  design  presented  here  is 
based  on  a  system  whose  continuous  state  and  output 
equations  can  be  defined  as: 

x(t)  =  Ax(t)  +  au  (t)  +  Dd(c)  (2-1) 

y  (t)  -  Cx  ( t )  +  Fu(t)  (2-2) 

where 

A  =  continuous-time  plant  matrix  ^nxn) 

B  =  cont inuous- t ime  control  matrix  (nxm) 

D  =  continuous-t ime  disturbance  matrix  (nxr) 

C  =  cont inuous- t ime  output  matrix  (pxn) 
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F  =  cont inuous- t ime  feed  forward  matrix  (pxm) 

The  dimensions  n,  m,  p  and  r  are  defined  as  the  number  of 
states,  inputs,  outputs  and  disturbances,  respectively. 
This  thesis,  however,  deals  entirely  with  the  case  where 
m«*p.  (i.e.  number  of  inputs  =  number  of  outputs). 


Professor  Porter's 

i 

techniques  that  produce 

iv 

Controller  which  can  be 

*  w > 

y 

k 

For  the  purpose  of  this 

the  digital  case  where  tne  system  state  and  output 
difference  equations  assume  the  form: 

x [ ( k+1 ) T ]  =  0x(kt)  +  4>u(kt)  4  Ad(kt) 
y(kt)  =  Tx(kT)  +  Yu(kT) 

where 

6  *  sampled-data  plant  matrix  (nxn) 

^  =  sampled-data  control  matrix  (nxm) 

A  *  sampled-data  disturbance  matrix  (nxr) 
r  3  sampled-data  output  matrix  (pxn) 

7  =  sampled-data  feed  forward  matrix 
T  =  sampling  period 

with 


(2-3) 

(2-4) 


9  =  e 

*■{ 


AT 


AT. . 
e  Bd  t 


(2-5) 

(2-6) 


A  -  ^eATDdt 

r  =  c 
y  =  f 

Figure  2-1  represents  the  Proportional  Plus  Integral 
Controller  for  the  digital  case.  The  discrete  integrator 
shown  in  this  figure  satisfies  Equation  (2-15).  For  this 
discrete  case,  samplers  are  assumed  to  be  placed  in  the 
feedback  loop  and  at  the  command  input  vector,  v.  A 
zero-order-hold  device  is  placed  after  the  second  summer  to 
produce  piece-wise  constant  inputs  to  the  plant.  An 
important  advantage  of  this  method  of  controller  design  is 
that  after  making  these  assumptions,  no  further  complex 
design  analysis  in  the  ’’Z"  domain  is  required. 


FIGURE  2-1 s  (PI)  Controller 
Using  this  type  of  system,  Professor  Porter  has 
developed  four  separate  design  procedures  with  which  to 
synthesize  a  Proportional  Plus  Integral  Digital  Controller. 


(2-7) 

(2-8) 

(2-9) 
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The  first  three  procedures  are  based  on  whether  the  plant 
involved  is  described  as  unknown,  known/regular,  or 
known/irregular.  The  fourth  procedure,  known  as  B*, 
provides  synthesis  of  the  digital  controller  to  provide 
enhanced  decoupling  and/or  elimination  of  any  undesirable 
initial  undershooting  of  the  outputs.  These  four  procedures 
are  briefly  described  in  the  following  paragraphs.  It 
should  be  noted  that  even  though  the  designs  discussed  are 
for  a  discrete  controller,  the  continuous  state  equation 
given  by  Equations  (2-1)  and  (2-2)  may  be  used  (Ref.  5). 

Unknown  Plants  (Ref .  6 ) 

In  many  cases,  very  little  is  known  about  the  dynamical 
properties  which  describe  the  state  equation  model.  In  this 
case  exact  representation  for  the  matricies  A,  B,  C,  D  and  F 
of  Equations  (2-3)  thru  (2-9)  are  not  available.  It  can 
therefore  be  very  expensive  and  time  consuming  to  develop 
such  a  model  when  one  does  not  already  exist.  At  the  same 
time,  the  state  transfer  function  matrix,  G(0),  can  be 
easily  determined  from  "off-line"  tests.  These  tests  must 
be  done  in  the  absence  of  the  distrubance  vector,  and  can  be 
accomplisned  only  if  the  open- loop  plant  is  asymptotically 
stable.  Also,  it  should  be  noted  that  when  the  state 
equations  are  available  and  the  designer  decides  to  use  this 
technique,  G(0)  can  be  calculated  by  Equation  (2-11).  In 
either  case  the  matrix  transfer  function  is 
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(2-10) 


G  (X)  =  C(XIn  -  A) _1B 


and 


U 


is 

•  i 


G(0) 


B 


(2-11) 


To  preserve  stability,  the  rank  of  G(0)  must  be  equal 
to  the  number  of  system  outputs.  It  is  further  required 
that  the  number  of  system  outputs  be  less  than  or  equal  to 
the  number  of  inputs  and  that  G(0)  have  no  transmission 
zeros  at  the  origin  or  in  the  right  half  of  the  s-plane. 

In  this  case,  the  control  law  equation  for  an 
error-actuated  Proportional  Plus  Integral  Digital  Controller 
assumes  the  form: 

u(kT)  »  T[€KQe(kT)  +€K1z(kT)]  (2-12) 

where 
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H-1 

(2-13) 

e(kT)  =  v(kT) 

-  y  ( kT) 

(2-14) 

z [ ( k+1 ) t]  =  z ( kT) 

+  Te (kT) 

(2-15) 

In  these  equations,  alpha  determines  the  ratio  of 
proportional  to  integral  control;  epsilon  is  a  scalar 
multiplier;  T  is  the  sampling  period;  v(kT)  is  the  command 
input  vector;  y(kT)  is  the  output  vector;  e(kT)  is  the  error 
vector;  z(kT)  is  the  integral  of  the  error  vector;  and  KQ 
and  Kj  are  given  by: 
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(2-16) 


Kq  =  orGT(0)  [G  (0)GT(0)  ]  _12 
Kx  =  GT(0) [G(0)GT(0) ]-12  (2-17) 

where  sigma  is  a  weighting  matrix  whose  scalar  diagonal 
elements,  <7.  ( t  =1 , 2 , 3 , . . . .  p)  are  chosen  by  the  designer.  If 
it  is  further  restricted  so  that  the  number  of  outputs 
equals  the  number  of  inputs,  as  is  the  case  for  this  thesis, 
then  G ( 0 )  is  a  square  matrix.  In  this  case  (G (0)G^(0) ]  1  = 

tGT(0)] 1  [G(O)]1  and  Equations  (2-16)  and  (2-17)  reduce 

to: 

Kq  =  aG(O)"1!  (2-18) 

Kx  =  G(0)_12  (2-19) 

Equations  (2-12)  thru  (2-19)  result  in  two  secs  of 
closed-loop  roots,  Z-^  and  Z2  : 

Z1  =  {A«C:|\In  -  In  -  TA  +  H  (T2 )  |  =  0}  (2-20) 

=  { Ae  C : | XI  -  I  +  T2  +  H ( T 3  >  |  =  0}  (2-21) 

where 

Z^  =  set  of  roots  due  to  plant  state* 

Z 2  =  set  of  roots  due  to  integrator  states 

I  =  nxn  identity  matrix 
n  J 

2 

HT  =  analytic  expression  whose  value  is  a 
function  of  t 

3 

HT  =  analytic  expression 
function  of 


whose  value  is  a 


Known  P  lant  (Re f .  5_)_ 

Kn  iwn  plants  fall  into  two  categories,  regular  and 
irregular.  For  a  plant  to  be  classified  as  known,  the 
matrices  A,  B,  C,  D,  and  F  of  Equations  (2-1)  and  (2-2)  must 
be  available.  To  design  the  controller  using  Porter's 
method  for  known  plants,  Equations  (2-1)  and  (2-2)  may  first 
be  put  into  the  form: 


V't) 

x2(t) 


A11  A12 
A21  A22 


and 


Y(t)  =  [C1  C2] 


xx  (t) 

Lx2(tij 


xL  (t) 

+ 

0 

■ 

_x2(t)_ 

.B2. 

r 

1 

u(t) 


(2-22) 


(2-23) 


where 


A11 

= 

A 

matr ix 

partition , 

(n-p)  x(n-p) 

A12 

= 

A 

matr ix 

partition , 

(n-p)  xp 

A21 

= 

A 

matrix 

partition, 

px  (n-p) 

A22 

= 

A 

matrix 

partition, 

pxp 

B2 

= 

non-zero 

B  matrix  partition,  fu 

C1 

= 

C 

matrix 

partition , 

px  (n-p) 

C2 

= 

C 

matrix 

partition, 

pxp 

If  Equations  (2-1)  and  (2-2)  cannot  be  put  into  the 
form  of  (2-22)  and  (2-23)  by  reordering  the  states,  they 
be  transformed  by  a  matrix  T  into  the  proper  form  by 


can 


choosing  T  in  x'(t)=Tx(t)  to  meet  this  requirement. 


T"1B  = 


-1 


0 

B, 

L  4J 


(2-24 ) 


Once  T  i3  found  that  satisfies  Equation  (2-24) ,  the  A 
matrix  of  Equation  (2-1)  can  be  transformed  using  the 
relationship 


A'  =  T-1AT 


and  the  C  matrix  is  transformed  using 


C'  *  CT 


Then  Equations  (2-1)  and  (2-2)  become: 


(2-25) 


(2-26) 


x' (t)  *  A'x(t)  + 

y  ( t)  =  c'x(t) 


o 

B. 


u(t) 


(2-27) 

(2-28) 


which  can  be  made  to  conform  to  Equations  (*.-22)  and  (2-23). 
It  should  be  noted,  however,  that  the  transformation  to  the 
form  of  Equations  (2-22)  and  (2-23)  is  not  necessary  in 
order  to  perform  the  design.  This  is  shown  in  the  example 
given  on  page  44  of  Reference  5.  However,  the  form  of 
Equations  (2-22)  and  (2-23)  are  used  in  the  discussion  to 
follow. 

If  the  original  system  contains  a  feedforward  matrix  F 
in  Equation  (2-2),  then  manipulation  of  the  original  system 


12 


so  as  to  remove  the  F  matrix  is  recommended.  This  allows 


the  system  to  be  represented  in  the  form  of  Equations  (2-22) 
and  (2-23).  This  can  be  accomplished  by  either  using 
outputs  which  do  not  contain  F  terms  (F  terms  are  caused  by 
accelerations  appearing  in  the  outputs)  or  by  augmenting  the 
original  system  with  the  servo  dynamics  yielding  new  state 
equations  which  do  not  contain  an  F  matrix. 

To  successfully  design  a  controller  using  Porter's 
techniques,  all  of  the  following  must  hold  true: 

(1)  The  pair  (A ,B)  must  be  a  controllable  pair 

(2)  The  pair  (A,C)  must  be  an  observable  pair 
A  B 

=  n+p 

-C  D 

(4)  p  _<  m  (for  an  irregular  design  p=m) 

In  addition  the  designer  must  check  the  rank  of  the 
first  Markov  parameter  to  determine  whether  the  plant  is 
regular  or  irregular.  The  first  Markov  parameter  is  defined 
as  CB=C  B  for  the  system  described  by  Equations 
(2-22)  and  (2-23).  If  the  system  is  not  in  the  form  of 
Equations  (2-22)  and  (2-23)  then  the  first  Markov  parameter 
is  defined  as  CB.  If  the  rank  of  this  parameter  is  equal  to 
the  number  of  outputs,  p,  then  the  plant  is  said  to  be 
"regular".  If  the  rank  is  less  than  p,  then  the  plant  is 


(3)  Rank 


Known/Re^ alar  P lants  (Ref.  5 ) 

If  the  first  Markov  parameter  has  full  rank  (i.e.  equal 
to  p)  the  plant  is  referred  to  as  being  "regular"  and  CB  is 
invertable.  When  this  is  true  the  system  can  be  represented 
by  the  block  diagram  of  Figure  2-2.  Although  a 
continuous-t ime  integrator  is  shown  in  Figures  2-2  and  2-3, 
a  discrete  integrator  of  the  form  shown  in  Figure  2-1  is 
intended . 


FIGURE  2-2:  Block  Diagram  of  Known/Regular  Plant/Controller 
Structure . 

For  this  system  the  control  law  equation  assumes  the 

form: 

u(kT)  =  (1/T) [K0e (kT)  +K1z(kT)]  (2-29) 

For  tracking  constant  commands  and  the  rejection  of  constant 
disturbance  the  controller  matrices  are  designed  such  that: 
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Kq  =  o[C0Bj_12 


2  2 


-1, 


' 1  "  fC2B2]  2  =  K0/a 


(2-30) 

(2-31) 


where 


S  =  diag  (  tf1  o2  . . . .  cp) 

<Tj  C  R 


(2-32) 

(2-33) 


The  error  vector  e(kT)  is  the  difference  between  the 
command  input  vector  v(kT)  and  the  output  vector  y(kT)  (see 
Figure  2-2).  Since  the  complete  system  is  augmemted  with  a 
vector  integrator,  the  steady-state  valua  of  the  error 
vector  is  zero  for  a  constant  command  vector  input,  and 
tracking  is  achieved  (Ref.  5). 

Furthermore,  this  design  theory  shows  that  increasingly 
'tight'  tracking  occurs  as  the  sampling  frequency,  f=l/T, 
approaches  infinity.  It  can  also  be  shown  (Ref.  5)  that  as 
T  approaches  zero  the  clcsed-loop  transfer  function  matrix 
G  ( X)  assumes  an  asymptotic  form: 


r  ( x )  =  r  ( x)  +  r  ( x) 


(2-34) 


where 


“lr 


r<M  =■  c0(Mn  -  in  -  taq )  tb0 


A 

r(x) 


-l. 


=  C2(XIp  ‘  Tp  +  9B2K0C2)  gB2K 


A„  = 


-ko~1ki 
A12C2  K0K1 


0 

A11  “  A2C2  C1 


(2-35) 

(2-36) 

(2-37) 
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B0  = 


A12C2 


-1. 


-1 


C0  =  tK0  Kx,  0] 


(2-38) 


(2-39) 


It  can  be  seen  from  Equations  (2-34),  (2-35)  and  (2-37) 
that  r(X)  contains  two  sets  of  roots,  associated  with  what 
are  called  the  "slow  modes".  These  two  sets  of  roots  are: 


“It 


ZL  =  {A  « C :  |  A  I  -  I  +  TKQ  J’K]_  |  =  0} 


(2-40) 


and 


-1, 


Z2  -  U*C:|AIn.p-In.p-TAu+TA12C2  "cj  =  0)  (2-4!) 

From  Equations  (2-34),  (2-36),  (2-38)  and  (2-39)  it  is  noted 

A 

that  the  poles  of  r(A)  contains  a  set  of  roots  associated 
with  the  "fast  modes".  This  set  of  roots  is: 


l3  =  { A « C : | X I  -  I  +  C2B2Kq|  =  0} 


(2-42) 


It  can  be  shown  (Ref.  5)  from  Equations  (2-37),  (2-38), 

and  (2-39)  that,  as  the  sampling  time  T  gets  smaller,  the 
slow  mode  roots  associated  with  the  set  of  roots 
become  uncontrollable,  and  the  slow  mode  roots  associated 
with  the  set  of  roots  Z2  (composed  of  the  transmission 
zeros)  become  unobservable.  Also,  the  fast  modes  remain 
both  controllable  and  observable.  Thus,  as  T  goes  to  zero, 
the  slow  modes  disappear  from  the  overall  transfer  function 
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T( X)  ,  and  only  the  fast  modes  remain: 

r(x)  =  r (\)  =  [m  -  ip  +  c2b2k0]-1c2b2k0  (2-43) 

The  r^ots  associated  with  the  fast  modes  are  chosen  by 
the  designer  to  lie  within  the  unit  circle  of  the  z-plane  by 
making 


=  ^  =  ^ i ag  ( <r A 1  (j„ •  •  •  on ) 


1 '  "2  P 


(2-44) 


and  is  found  using  Equation  (2-30).  Combining 
Equations  (2-43)  and  (2-44)  yields: 


r(x)  = 


=  diag 


1  p  +  £i 


X  -l+<r 


X-1+<T„ 


X-l+<r, 


(2-45) 


and  thus  Equation  (2-45)  shows  that  decoupling  is  achieved 
in  the  asymptotic  case. 

Transmission  zeros  of  the  system  are  a  suoset  of  the 
slow  modes  and  can  be  found  using  Equation  (2-41)  if  the 
original  system  has  the  form  of  Equations  (2-22)  and  (2-23). 
If  this  is  not  the  case,  the  transmission  zeros  can  be  found 
using  a  computer  program  described  in  Reference  7. 
Transmission  zeros  can  be  a  problem  since  Porter's  method 
requires  them  to  lie  within  the  unit  circle  in  the  z-plane. 
If  they  are  located  on  or  outside  the  unit  circle  then  a 
stable  controller  maynot  be  achievable. 
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Ktiown/I  r  regular  P  lants  (Ref .  5 ) 

If  the  first  Markov  parameter  is  rank  deficient  (i.e. 
not  equal  to  p)  then  the  plant  is  said  to  be  "irregular". 

In  this  case  the  product  CB  is  not  invertable  and  it  becomes 
necessary  to  introduce  measurement  equations  to  compensate 
for  this  rank  deficiency.  When  this  happens  a  new  feedback 
vector  is  defined  as: 


w  (t.) 


tC1  +  MA1X 


c2  +  ma121 


X1  (fc) 

x2(t) 


f2] 


XL  (t) 
X2(t) 


(2-46) 


(2-47) 


where 


M  =  Measurement  Matrix,  px(n-p) 

The  matrix  'M*  is  a  measurement  matrix  which  is  chosen  such 
that  the  matrix  product  FB  has  full  rank.  With  the  addition 
of  the  measurement  matrix  the  system  can  now  be  represented 
by  the  Block  Diagram  of  Figure  2-3. 


FIGURE  2-3:  Block  Diagram  of  Known/Irregular  Plant/Controller 
S  tructure . 

The  measurement  matrix  must  be  chosen  by  the  designer. 
Only  through  proper  choice  of  this  matrix  can  decoupling  of 
the  outputs  in  the  asymptotic  case  occur.  Reference  8 
discusses  a  method  for  choosing  a  measurement  matrix  which 
allows  for  decoupling.  The  method  discussed  in  this  paper 
works  only  if  the  system  can  be  put  in  the  form  of  Equations 

(2-22)  and  (2-23).  Once  in  this  form  a  B*  matrix  is  formed 

from  the  following  equation: 

C  TA  d  A  1 

1  >  12  (2-48) 

CmTAlld,”fl12. 

where  m  is  the  number  of  control  inputs,  C-T  is  the  ith  row 

of  C1  and  dj  =min  ( j  :C.T  A^  A12”0'  3  =  0  • 1  • 2  •  •  •  n-1 )  or  dj  =n-l 

T  j 

if  C*  A^  Ai2=^  f°r  Next  F2  formed  using: 

F2  =  C2  +  MAI2  (2-49) 


and  by  comparison  of  Equation  (2-48)  to  (2-49)  the  m 
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elements  ace  chosen  based  upon  the  non  zero  elements  of  B*. 
It  should  be  noted  that  additional  constraints  are  put  on 
the  elemer's  of  the  M  matrix  depending  on  the  rank  of  B*  and 
the  row  dimension  of  as  compared  to  row  dimension  of 
and  column  dimension  of  A^.  For  a  complete  understanding 
of  the  choice  of  the  M  matrix  elements  the  reader  must  refer 
to  Reference  8.  Also  it  should  be  pointed  out  that  proper 
choice  of  the  measurement  matrix  elements  is  critical 
because  it  affects  the  location  of  the  transmission  zeros  of 
the  system.  These  transmission  zeros  are  now  determined  by: 

ZT  »  U^lAI^-I^p-A^A^F-f^!  =  0)  (2-50) 

and  as  with  regular  plants,  the  transmission  zeros  must  lie 
in  the  left  half  of  the  s-plane  or  within  the  unit  circle  in 
the  z-plane.  Once  again,  if  the  system  is  not  represented 
as  in  Equations  (2-22)  and  (2-23),  then  Equation  (2-50)  can 
not  be  used  to  determine  the  transmission  zeros. 

Once  the  measurement  matrix  M  is  chosen  and  F  and  F 
are  calculated,  the  control  law  is  governed  by  the  form: 


u(kt)  -  (1/T) [K0e (kT)  +  K1z(kT)] 


(2-51) 


where 


Kq  =  «(F2B2)’12 

K1  -  «f2b2)_1s 


(2-52) 

(2-53) 


Note  that  C-,  in  Equations  (2-30)  and  (2-31)  is  replaced  by  F 
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in  Equations  (2-52)  and  (2-53),  and  that  does  have  full 

rank,  by  proper  choice  of  M,  and  is  therefore  invertable. 

Figure  2-3  shows  that  the  error  vector  e(kT)  is  now  given 
by: 


e ( kT)  =  v(kT)  -  w(kT)  (2-54) 

It  can  also  be  shown,  as  before,  that  for  constant  inputs, 
the  steady-state  value  of  w(kT)  equals  y(kT),  so  that  the 
system  still  tracks  the  constant  command  input,  as  desired. 

As  with  Known/Regular  plants,  it  can  be  shown  that  as 
the  sampling  time  T  approaches  zero  the  asymptotic 
closed-loop  transfer  function  G(X)  assumes  the  form: 

T(X)  =  F( X)  +  ?(X)  (2-55) 

where 


f  (X) 

co 

1 

C 

M 

✓< 

Xn  '  TA0 

)  TB0 

(2-56) 

A,  . 

<xip- 

_  1 

r(xj 

co 

h 

?  -  a4> 

~B2K0 

(2-57) 

'-Ko'1 

Ki 

0 

Ao 

-l 

-1 

(2-58) 

A12F2 

K0K1 

A11  A12F2  F1 

- 

0 

* 

B  9 

— 

-1 

(2-59) 

.A12F2 

[c 

-1 

-1 

-] 

C0 

2K2  K 

0 

K1  ' 

^1  -  C2F2 

(2-60) 

A4 

= 

B2K0F2 

(2-61) 

21 


Ac  expected,  these  equations  are  simiular  to  Equations 
(2-34)  thru  (2-39)  which  were  developed  for  Known/Regular 
plants.  Because  of  this  singularity  it  can  be  seen  that  the 
two  sets  of  "slow  mode"  roots  are  now: 

Zx  =  {X«C: |Xlp-Ip+TK0"1K1|  =  0}  (2-62) 

and 


z2  -  {X«C:|XIn.p-In.p-TA11+TA12F2-  FjJ  -  0)  (2-63) 


and  that  the  "fast  mode"  roots  become: 


Z 3  =  {XeC:|XI  -I  +F2B2K0|  =  0}  (2-64) 

Again,  analogous  to  the  Known/Regular  plant,  it  can  be 
shown  (Ref.  5)  that  as  the  sampling  time  T  gets  smaller  the 
slow  modes  associated  with  the  Z1  roots  become 
uncontrollable  and  thus  have  no  contribution  to  the  transfer 
function.  But,  unlike  the  case  of  the  Known/Regular  plant, 
the  slow  modes  associated  with  the  Z2  roots  remain  both 
observable  and  controllable  due  to  the  extra  measurements 
generated  by  the  M  matrix.  The  fast  modes  associated  with 
roots  (Zg) ,  again  remain  both  controllable  and  observable  as 
T  gets  smaller. 

Hence,  it  is  evident  that  as  f  approaches  infinity  the 
transfer  function  matrix  G(X)  of  the  d iscre te-t ime  tracking 
system  assumes  the  asymptotic  form: 


-1 
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r(X)  =  (O-CjFj  Fp  IXIn-p'In-p"TAll+TA12F2  F1I  TA12F2 

+C2F2~1[XVIp+F2B2K01"lF232K0  <2-65> 


The  first  terra  of  Equation  (2-65)  contains  the  slow  mode 
roots  (transmission  zeros)  anu  the  second  term  contains  the 
fast  mode  roots. 

The  fast  mode  roots  are  chosen  by  the  designer,  using  a 
method  simiular  to  that  used  for  Known/Regular  plants: 


^2® 2^0  =  ^2  +  ^12®2^ 0  =  d i ag  (  or^  ,  •  •  (Jp)  (2—66) 

These  roots  must  lie  within  the  unit  circle  in  the  z-plane. 
The  transmission  zero  locations  are  also  set  by  the  designer 
in  Equation  (2-63)  and  are  chosen  by  selection  of  the 
measurement  matrix  M. 

It  should  be  noted  that  the  transfer  function  matrix 
given  in  Equation  (2-65)  does  not  always  result  in  decoupled 
outputs  and  that  the  measurement  matrix  M  plays  a  large  role 
in  determining  decoupling.  Again  the  reader  is  referred  to 
Reference  8  for  help  in  picking  a  measurement  matrix  that 
allows  decoupling  of  the  outputs. 
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B*  Plants  (Ref.  9) 

For  plants  in  which  the  three  previous  designs  do  not 
satisfactorily  decouple  the  output  and/or  undesirable 
initial  undershooting  occurs  in  the  time  response,  this 
fourth  design  method  is  available.  A  new  matrix,  B*  is 
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determined  using 


B*  = 


rClTAdlB 

=2V2B 


c  TAd?B 
P 


(2-67) 


where  m  is  the  number  of  control  inputs,  Cp  is  the  ich  row 
of  C,  and  dp  -  min ( j :Cp^AiB=0 , j=0 ,1 , . . . n-i )  or  dp  =  n-1  if 
Cp^A*B  =0  for  all  j.  Note,  however,  that  if  B*  is  rank 
deficient  (i.e.  B*  is  singular)  this  design  method  cannot  be 
used. 

When  B*  has  full  rank,  the  controller  matrices  are 
determined  by: 


K0  =  (B*)  121 
Kx  =  [G  (0 )  ]"12, 


(2-68) 

(2-69) 


where 21  and  22  are  both  selected  as  diagonal  (pxp) 
matrices.  Once  Kq  and  have  been  determined  by  Equations 
(2-66)  and  (2-67) ,  the  control  law  is  governed  by  the 
equation: 


u  (kt)  =  (1/T)  [KQe  (kT)  +  K]_z(kT)) 


(2-70) 


with 


e (kT)  =  v(kT)  -  Y (kT) 


(2-71) 


Analysis  of  this  method  reveals  that  Kn  shapes  the 
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initial  output  response  (i.e.  no  undershoot)  while  K1  shapes 
the  steady-state  outputs  (i.e.  no  offset).  It  is  also 
important  to  note  that  G(Q)  and  B*  can  both  be  obtained  from 
1  1  '  tests  of  the  open— loop  plant  by  using  the 

unit-step  response  matrix 

S  (T)  =  jfceAtBdt  (2-72) 

such  that 

S (®)  =  G(0),  as  in  the  unknown  case  (2-73) 

S (0 )  =  B*  (2-74 ) 


Summary 

This  chapter  presents  a  summary  of  the  four  digital 
controller  design  methods  developed  by  Professor  Brain 
Porter.  The  highly  interactive  computer  program  MULTI 
(Ref.  4)  implements  the  four  design  methods  discussed.  It 
should  also  be  pointed  out  that,  although  this  chapter 
emphasises  the  digital  approach,  analog  controllers  can  also 
be  designed  using  the  same  overall  concepts.  Chapter  III, 
the  next  chapter,  gives  a  description  of  the  KC-135A 
aircr  ft  along  with  a  discussion  of  the  equations  of  motion, 
sign  convention,  and  axis  system  used. 


Ill . 


KC-135A  AIRCRAFT  MODEL 


Introduction 

This  thesis  deals  with  the  KC-135A  aircraft.  A  linear 
six  degree-of-f reedom  model  is  generated  for  each  of  the 
three  flight  conditions  chosen.  This  results  in  a  set  of 
three  linearized  constant-coefficient  equations  of  motion, 
each  one  representing  the  aircraft  operating  within  small 
perturbations  about  its  respective  equilibrium.  For  the 
purpose  of  this  thesis  the  three  flight  conditions  are 
selected  so  as  to  represent  the  aircraft  over  a  wide  range 
of  flight  operating  conditions.  This  has  been  done  so  that 
one  controller,  however  chosen,  could  be  tested  for 
robustness  at  the  three  distinctly  different  flight 
conditions.  Thus  the  flight  conditions  chosen  include  a 
high  altitude  high  speed  cruise  of  MACH  0.77  at  45,000  feet 
a  medium  altitude,  heavy  weight  cruise  of  MACH  0.77  at 
28,500  feet;  and  a  landing  condition  of  MACH  0.21  at  sea 
level . 

This  chapter  provides  a  basic  description  of  the 
KC-135A  aircraft  and  overviews  a  general  discussion  of  the 
linearized  equations  of  motion,  sign  convention,  and  axis 
system  used  in  this  thesis. 

Aircraft  Description  (Ref.  10 ) 

The  KC-135A  is  a  four  engine  jet-powered  tanker/cargo 
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airplane.  The  swept  wing  is  mounted  low  on  the  fuselage  at 
an  inc idence  of  2  degrees  and  is  tailored  for  high  subson ic 
cruise  speeds.  The  aircraft  has  a  basic  weight  of 
approximately  106,000  pounds,  depending  on  equipment 
installed,  and  a  maximum  gross  weight  of  287,000  pounds. 

All  control  surfaces,  except  the  spoilers,  are 
aerodynamically  balanced  and  operated  by  means  of  control 
tabs.  A  hydraulically  boosted  rudder  is  installed  on  all 
aircraft.  The  lateral  control  system  is  composed  of 
integrated  aileron  and  spoiler  control  surfaces.  The 
spoilers  may  also  be  used  as  speed  brakes  when  operated 
symmetrically.  Movement  of  the  inboard  ailerons  causes  a 
corresponding  movement  of  the  outboard  ailerons  if  the  wing 
flaps  are  extended  beyond  the  23  degree  point.  Such  is  the 
case  for  the  flight  condition  representing  the  landing  phase 
in  this  thesis.  If  the  wing  flaps  are  up,  a  lockout 
mechanism  prevents  the  outboard  ailerons  from  moving. 

Lateral  trim  of  both  the  rudder  and  ailerons  is  accomplished 
manually  by  rotating  a  trim  wheel  which  positions  trim  tabs 
on  respective  control  surfaces. 

Longitudinal  control  is  provided  by  an  all  moveable 
horizontal  stabilizer  and  elevator  system.  The  stabilizer 
position  is  set  by  a  trim  wheel  which  can  be  operated 
electrically  or  manually.  All  three  flight  conditions  used 
in  this  thesis  assume  a  horizontal  stabilizer  setting,  which 
results  in  no  elevator  deflection  required  to  maintain  that 


flight  condition.  It  should  also  be  noted  that  no  flight 
control  surface,  either  lateral  or  longitudinal,  is  modified 
for  the  purpose  of  this  thesis  and  all  function  as  described 
by  the  current  Tech.  Order  specification  at  the  time  of  this 
thesis  (Ref.  11) 

The  four  Pratt  and  Whitney  J57-P-59W  or  -43WB  engines 
are  mounted  individually  below  the  wing  on  forward  swept 
struts.  The  engines  are  each  rated  at  12845  pounds  thrust 
for  a  standard  day  at  sea  level  (15  deg.  C,  29.29  inches  of 
mercury).  Other  relative  geometric  data  as  found  in 
Reference  10  is  as  follows: 


Characteristic 

Symbol 

Dimension 

Fuselage 

F 

128.83  ft 

Wing  Area 

S 

2433  sq.ft 

Wing  Spa” 

b 

130.83  ft. 

Wing  M..  •'  . 

c 

20.16  ft. 

Distance  from  25% 
Wing  M.A.C.  to 

25%  Horizontal 

Tail  M.A.C. 

/T 

61.39  ft. 

The  abbreviation,  M.A.C  ,  s’  ds  for  mean  aerodynamic  chord. 
Further  data  for  the  KC-135A  aircraft  can  be  found  in 
References  10  and  11. 

System  Models 

The  aircraft  models  in  this  thesis  are  developed  using 
six  degree-of-f reedom  equations  of  motion.  The  equations 
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are  assumed  to  have  no  coupling  between  motion  in  the 
lateral  directional  plane  and  the  longitudinal  plane. 
Therefore,  the  original  six  differential  motion  equations 
can  be  decoupled  into  two  sets  of  three  equations  each.  One 
set  describes  motion  of  the  aircraft  in  the  longitudinal 
plane  or  plane  of  symmetry  while  the  other  describes  motion 
only  in  the  lateral  and  directional  plane  or  motion  out  of 
the  plane  of  symmetry.  Appendix  A  presents  the  development 
of  the  equations  of  motion  in  state  form  along  with  othe. 
assumptions  made  about  them. 

Control  inputs  available  include  two  lateral  controls, 
the  rudder  (j  )  and  the  ailerons  and  spoilers  ({  ).  6  is 

modeled  as  the  control  wheel  movement  and  its  limits  are  set 
at  +  90  degrees.  A  control  wheel  movement  would  correspond 
to  a  combined  aileron  and  spoiler  displacement.  Movement  of 
the  control  wheel  within  its  limits  would  result  in  both 
aileron  and  spoiler  movements  within  their  limits.  is 

rudder  displacement  and  its  limits  are  set  at  +  17  degrees. 
Longitudinal  controls  include  the  elevator  (6e )  ,  the  speed 
brakes  (6^ )  ,  and  the  thrust  (i^).  SQ  is  elevator 
displacement  and  its  limits  are  set  at  +  25  degrees.  It 
should  be  noted  that  any  constant  elevator  displacement 
represents  an  out  of  trim  condition  and  would  normally  be 
trimmed  off  by  an  appropriate  horizontal  stabilizer 
position.  However,  this  is  not  possible  with  the  model 
representation  used  in  this  thesis.  <5^  is  speed  brake 
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deflection  which  is  symmetrical  spoiler  deployment  and  has  a 
limit  of  0  to  60  degrees.  5^.  is  modeled  in  terms  of  100%  of 
available  thrust  and  limits  vary  depending  upon  altitude  and 
gross  weight. 

Since  the  equations  are  decoupled  and  the  theory 
requires  that  the  number  of  inputs  equals  the  number  of 
outputs,  two  lateral  and  up  to  three  longitudinal  outputs 
are  chosen.  The  lateral  outputs  are  the  roll  angle  <p  and 
the  sideslip  angle  0.  The  longitudinal  outputs  depend  on 
the  maneuver  being  performed  and  the  controls  available. 

This  is  necessary  since,  in  the  case  of  a  constant  airspeed 
climb,  the  speed  brakes  are  not  used  and  are  removed  from 
the  model.  Chapter  V  gives  a  more  detailed  description  of 
the  inputs  and  outputs  used  for  each  maneuver. 

S ign  Convention  and  Axis  System 

The  sign  conventions  for  forces  and  moments,  as  used 
throughout  this  thesis,  are  shown  in  Figures  3-1,  3-2,  and 
3-3.  Figure  3-2  shows  lateral  sign  conventions  while  Figure 
3-3  shows  longitudinal  sign  conventions. 

Rudder  (5f )  :  Rudder  deflection  to  the  left  is  defined 
as  positive.  This  produces  a  positive  0  ,  positive  V, 
negative  N,  and  negative  R. 

Control  Wheel  (5..,)  :  Control  wheel  deflections  to  the 

. .  “'“■■■  w 

right,  which  cause  right  aileron  up  and  right  spoiler  up 
along  with  left  aileron  down  is  defined  as  positive.  This 
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P  -  Roll  Angular  Velocity 

Q  -  Pitch  Angular  Velocity 

R  -  Yaw  Angular  Velocity 

X,Y,Z,  -  Aerodynamic  Force 
Components 

ufv,w  -  Velocity  Components 


-  Rolling  Moment 

-  Pitching  Moment 

-  Yawing  Moment 

-  Roll  Angle 

-  Pitch  Angle 

-  Yaw  Angle 


FIGURE  3-1:  Sign  Convention 
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produces  a  positive  L,  positive  0,  and  positive  P.  Note, 
that  this  definition  of  positive  aileron  is  not  standard, 
but  does  conform  with  Reference  12. 

Elevator  (6^):  Control  column  movement  forward  which 
caused  down  elevator  is  defined  as  positive.  This  positive 
elevator  produces  a  negative  $  ,  negative  M,  and  negative  Q 

Speed  Brakes  (5^)  :  Spoilers  when  used  symmetrically 
are  defined  as  speed  brakes,  which  when  deflected  positive 
are  up. 

Thrust  (6)1  Thrust  is  modeled  as  %  of  available 
thrust  and  a  positive  6^  calls  for  an  increase  in  thrust. 

Body  Axis  System  This  thesis  uses  the  body  axis 
system  shown  in  Figure  1.  It  should  be  noted  that  Appendix 
A  takes  the  nondimensional  stability  axis  derivatives  found 
in  Reference  13  and  converts  them  to  dimensional  body  axis 
values. 


Summary 


This  chapter  presents  a  physical  description  of  the 
KC-135A  aircraft.  The  system  models,  as  developed  in 
Appendix  A,  are  overviewed.  Also,  sign  conventions  and  the 


axis  system  used  in  this  thesis  are  defined.  Chapter  IV, 
the  next  chapter,  presents  a  computational  time  delay 
modification  made  to  the  computer  program  MULTI  by  this 


author . 
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IV.  TIME  DELAY  MODIFICATION 
to~USmpUter  PROGRAM  MULTI 

Introduction 

The  core  of  this  thesis  is  developed  around  the  use  of 
a  computer  program  MULTI  (Ref.  4)  which  was  developed  to 
design  and  simulate  the  techniques  of  Professor  Brain  Porter 
discussed  in  Chapter  II.  MULTI  is  a  highly  interactive 
program  designed  to  give  the  user  the  opportunity  to  both 
design  a  multivarable  controller  and  simulate  it  with  the 
help  of  a  computer.  The  computer  provides  a  means  for 
multiple  iterations  of  the  design  and  simulation  to  allow 
the  user  to  achieve  a  suitable  design.  For  an  understanding 
of  the  program's  development  and  structure,  as  well  as  the 
basic  User  and  Programmer  Manual,  the  reader  is  referred  to 
Reference  4.  This  reference  provides  the  reader  with  a 
basic  setup  of  the  program.  However,  MULTI  has  gone  through 
several  revisions  and  to  find  information  on  the  most 
current  status  of  MULTI  the  reader  should  contact  Professor 
John  J.  D'Azzo,  Deputy  Department  Head,  Department  of 
Electrical  Engineering,  Air  Force  Institute  of  Technology, 

Wr  ight-Patterson  Air  Force  Base,  Ohio,  45433. 

The  remainder  of  this  chapter  presents  a  discussion  of 
the  time  delay  modification  made  to  MULTI  by  this  author. 

Theory  of  Computational  T ime  Delay  (Ref .  14 ) 

The  simulation  portion  of  the  computer  program  MULTI 
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provides  the  designer  with  a  method  of  determining  how  well 
his  design  will  perform.  The  simulation,  up  until  the  time 

of  this  author's  modification,  made  the  assumption  that 
there  is  no  computational  time  delay  for  generation  of  the 
inputs  to  the  plant.  In  some  cases  this  assumption  may  not 
be  valid.  Therefore  it  became  necessary  to  provide  the  user 
an  option  which  includes  a  computational  delay  in  the 
simulation  if  desired.  At  the  same  time  that  the 
computational  time  delay  was  added  to  MULTI,  a  compensator 
of  the  form  presented  in  Reference  14  was  also  added. 

Indeed,  if  no  such  compensation  is  provided  the  resulting 
tracking  system  may  have  either  very  poor  performance  or 
even  be  unstable. 

The  controller  of  Chapter  II  generates  a  signal  of  the 

form: 

s(kT)  =  f [KQe (kT)  +  KLz(kT)]  (4-1) 

where  f  is  the  sampling  frequency,  e(kT)  is  the  error 
vector,  z(kT)  is  the  integral  of  the  error  vector,  and  Kq 
and  K1  are  controller  matrices  of  the  form  discussed  in 
chapter  II.  Reference  14  proposes  a  computational  time 
delay  equation  of  the  form: 

r ( kT)  =  s (kT)  -  r [ ( k-1 ) T]  (4-2) 

Thus,  with  a  computational  time  delay  of  one  sampling 
period,  the  digital  controller  is  required  to  generate  the 
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input  vector 


u (kT)  =  r [  (k-l)T] 


(4-3) 


For  computational  time  delays  of  more  than  one  sampling 
period  Equation  (4-2)  becomes 


r (kT)  =  s(kT)  -  r (  (k-1) T]  -  r [ ( k-2 ) T ] 


_ r ( (k-x) T]  (4-4) 


and  the  new  input  vector  is  now  given  by 


u(kT)  »  r  [  (k-x) T] 


(4-5) 


where  x  is  the  number  of  sampling  periods  equal  to  the 
computational  delay.  Figure  4-1  gives  a  generalized  block 
diagram  for  this  computational  time  delay  and  compensation. 


i(kTW  r(kT) 


r((lc-l)T) 


r((k-i)T) 


r((k-x)T) 


FIGURE  4-1:  Computational  Time  Delay  VJith  Compensation 


Reference  14  develops  the  theory  for  the  new  system 
with  this  computational  time  delay  and  the  compensation 
represented  by  Equations  (4-3)  and  (4-2)  respectively.  This 
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theory  shows  that  the  time  delay  has  no  affect  on  the  "slow" 
modes  as  f  approaches  infinity  but  the  "fast"  modes  are 

affected.  In  fact,  it  can  be  shown,  Reference  14,  that  as  f 
approaches  infinity,  if  no  time  delay  compensation  is  used, 
at  least  some  of  the  "fast"  modes  become  unstable.  For  the 
complete  analysis  of  this  system  the  reader  is  referred  to 
Reference  14. 


Summary  of  MULTI  Code  Changes 

To  accommodate  this  time  delay  a  few  simple 
modifications  were  made  to  the  existing  Fortran  code  for 
MULTI.  First,  under  OPTION  #23  the  user  is  given  the 
opportunity  to  set  the  amount  of  computational  delay 
desired.  Computational  delay  is  limited  to  an  integer 
multiple  from  1  to  5  of  the  sampling  time  used.  This 
integer  is  passed  throughout  the  program  as  "NCD"  and  is  now 
entered  along  with  the  sampling  time  in  OPTION  #23.  The 
program  still  allows  for  no  computational  delay  if  NCD  is 
set  to  zero.  The  code  under  OPTION  23  reads  as  follows: 


2ii=C  *»*»t«»*»»*t*»*H****»***»tm***ttt»»t»*t*mttmttt»tt»Ht** 

218=C  OPTION  123 

228=  2123  PRINT*,  'THIS  OPTION  SETS  THE  SAMPL7N6  TIME  AND’ 

238*  PRINT’ (A/)’,  ’  COMPUTATIONAL  DELAY  TIME’ 

248*  PRINT*, ’ENTER  THE  SAHPLIN8  TIME  >’ 

258*  READ*,SAMPT 

268*  PRINT*,  ’  ’ 

278*  PRINT*,  'COMPUTATIONAL  DELAY  IS  AN  INTE6ER  MULTIPLE’ 

288*  PRINT*,  ’FROM  I  TO  5  OF  THE  SAMPLING  TIME.  USE  ZERO  IF’ 

290*  PRINT’ (A/)’,  ’  NO  DcLAY  IS  DESIRED' 

300=  1851  PRINT*,  'ENTER  THAI  INTEGER  >’ 


310=  READ*, NCD 

320=  IF  (NCD.LE.5)  SO  TO  1852 

330=  PRINT#,  ’  ’ 

340*  PRINT#,  ’COMPUTATIONAL  DELAY  MUST  BE  AN  INTEGER  MULTIPLE’ 

331*  PRINT#,  ’FORM  1  TO  5  OF  THE  SAHPLIN6  TIME.  USE  ZERO  IF' 

360=  PRINT’ (A/)',  ’  NO  DELAY  IS  DESIRED’ 

370=  60  TO  1631 

380*  1852  IFLA8(23)*1 
390=  GO  TO  8017 


When  OPTION  23  is  executed  lines  220  thru  250  prompt  the 
user  with  the  following  statements: 


THIS  OPTION  SETS  THE  SAMPLING  TIME  AND 
COMPUTATIONAL  DELAY  TIME. 

ENTER  THE  SAMPLING  TIME> 


Here  SAMPT  is  entered  and  the  program  continues  with  lines 
270  thru  300  and  the  user  is  prompted  with: 


COMPUTATIONAL  DELAY  IS  AN  INTEGER  MULTIPLE 
FROM  1  TO  5  OF  THE  SAMPLING  TIME.  USE  ZERO  IF 
NO  DELAY  IS  DESIRED. 

ENTER  THAT  INTEGER> 


The  integer  entered  here  is  stored  as  NCD  and  is  passed  to 
various  portions  of  the  program  by  the  common  statement: 


COMMON  /B  12B/  ST ,TT , SAMPT , NCD 


Lines  320  thru  370  of  the  code  under  OPTION  #23  prevent  the 
user  from  entering  a  value  for  NCD  which  is  out  of  the  range 
from  0  to  5.  This  is  necessary  to  prevent  overflow  of  the 
set  aside  storage  space  under  OPTION  #26,  which  is  explained 
later . 
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Code  also  was  also  added  to  OPTION  #29  to  allow  NCD  to 
be  read  in  from  a  memory  file  and  to  OVERLAY  (14,0)  to  write 
NCD  to  a  memory  file  when  OPTION  #99  is  executed.  The  read 
statement; 


RE  AO  v  NCD 


was  added  to  OPTION  #29  and  the  write  statement; 


WRITE  (DAT3 , *)  NCD 


was  added  to  OVERLAY  (14,0).  Under  OPTION  #123  TO  #126  the 
code  was  changed  to  read  as  follows. 


400-C  iHHHHHHiHHHHHHHHHHHHHHHHHHHHHHHtit 
4U*C  OPTIONS  1123  TO  1123 

421*  2123  PRINT  * (A, 191 , Fli. 5) ’ ,  '  SANPLINS  TIME. . . ' , SAHPT 
431*  CD=NCD*SANPT 

44#=  PRINT’ (A, 13X, F1#.S) ’ ,  ’  COMPUTATIONAL  DELAY...’, CD 
43S*  PRINT  ’ (A, 172,F1«.5)J ,  ’  SIMULATION  TIHE...’,TT 
441*  PRINT  ’ (A,11X,F1#.5) ’  CALCULATION  STEP  SIZE...’, ST 
471*  60  TO  8f IS 


Line  430  computes  the  computational  delay  time  using 
the  equation  CD  =  (NCD ) * (SAMPT) .  Line  440  adds  the  print 
out  of  CD  to  this  option  so  that  the  user  sees  the  following 
example  print  out  for  one  delay  period  and  T=  0.050  seconds 
whenever  OPTION  #123  TO  #126  is  called: 


SAMPLING  TIME  . . .  .  GdO 
COMPUTATIONAL  DELAY  ...  .050 
SIMULATION  TIME  ...  20.000 
CALCULATION  STEP  SIZE  ...  .050 


When  the  simulation  OPTION  #26  is  executed  several 
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additional  steps  occur  it  the  user  has  selected  a 
computational  delay  under  OPTION  #23.  First  a  storage 


matrix  UD  is  dimensioned  as  10  by  6.  This  matrix 
accommodates  a  delay  of  up  to  5  sampling  periods,  thus  the 
restriction  set  in  OPTION  #23  on  the  size  of  NCD.  The 
matrix  UD  is  then  initialized  to  zero  and  the  user  is 
reminded  that  a  computational  delay  is  being  included  in  the 
s imulat ion . 

The  code  that  accounts  for  the  computational  delay  is 
added  immediately  following  the  formation  of  the  control  law 
in  the  existing  program.  As  mentioned  earlier,  this  code 
can  accommodate  up  to  5  periods  of  delay  and  provide  for  the 
compensation  given  by  Equation  (4-4) .  This  code  as  it  now 
appears  in  MULTI  is  listed  below  and  resides  in  the  main 
loop  which  is  executed  for  each  sampling  period. 


S2I-C - DELAY  OF  U  BY  NCD  TINES - 

531*  IF  (NCD.EQ.I)  GO  TO  12U 
541*  DO  1211  J*1,MCD 

351*  DO  1212  1*1, H 

5GB*  U(I)*U(I)-UD(I, J) 

371*  1212  CONTINUE 
581*  1211  CONTINUE 
391*  DO  1216  1*1, H 
611*  1216  UD(I,NCD+1)*U(I) 

611*  DO  1207  1*1, N 
620*  1207  U(I)*UD(I,1) 

630*  DO  1208  J*t,NCD 

640*  DO  1209  1*1, H 

630*  UD(l,J)*UD(!,J+l) 

660*  1209  CONTINUE 
670*  120B  CONTINUE 

Line  530  is  an  IF  statement  that  checks  to  see  if  the 
integer  NCD  equals  zero,  which  would  represent  a 
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computational  delay  of  zero.  If  this  is  true  the  program 
skips  the  next  fourteen  lines  of  code  and  no  delay  is 

encountered.  Investigation  of  the  remaining  list  of  code 
reveals  the  following  fact  about  the  matrix  UD.  This  matrix 
is  used  to  store  the  input  vectors  that  are  generated  by 
Equation  (4-2)  and  thus  acts  as  th-.  mechanism.  The 

earest  way  to  understand  this  code  :  co  take  a  simple 
example  where  NCD  equals  one  and  to  walk  through  the 
process.  The  first  time  through  the  main  loop  a  vector 
s(kT)  is  formed  as  given  by  Equation  (4-1).  Since  NCD 
equals  one,  lines  540  thru  580  subtract  from  s(kT)  the 
vector  that  resides  in  the  first  column  of  UD  to  form  r(kt) 
in  accordance  with  Equation  (4-2).  The  first  column  of  UD 
always  contains  the  vector  r[(k-l)T],  Next,  lines  590  and 
600  store  r(kT),  given  by  Equation  (4-2),  in  the  NCD+1 
column  of  UD.  Lines  610  and  620  assign  the  input  vector 
u(kT)  from  the  first  column  of  UD  which  corresponds  to 
performing  Equation  (4-3)  and  thus  a  time  delay  of  one 
sampling  period  is  encountered.  Finally,  lines  620  thru  670 
move  the  second  column  of  UD  into  the  first  column  so  that 
the  program  is  ready  for  the  next  time  the  main  loop  is 
executed.  It  is  important  to  note  that  the  computer  code 
does  not  use  the  symbols  s(.)  and  r(.)  in  order  to  save 
storage  space. 

If  NCD  equals  3  it  can  be  shown  that  the  matrix  UD  is 
composed  of  the  following: 


UD  = 


r[(k-3)T]  ri(k-2)T]  r[(k-l)T]  r[kT] 


from  this  it  is  easy  to  see  how  UD  is  used  in  execution  of 
Equations  (4-4)  and  (4-5)  and  that  UD  indeed  acts  as  the 
delay  mechanism. 


Summary 

Fortran  code  changes,  as  discussed  in  this  chapter, 
were  made  to  the  computer  program  MULTI  to  incorporate  a 
computational  time  delay.  Debugging  of  the  code  became 
necessary  because  unexpected  results  were  initially  obtained 
when  example  problems  were  worked.  An  error  was  found  in 
the  existing  code  and,  once  corrected,  favorable  results  are 
obtained  both  with  and  without  computational  time  delay. 


V.  MULTIVARIABLE  CONTROL  LAW  DESIGN 


Introduction 

This  chapter  presents  the  procedures  used  in  the 
development  of  multivariable  tracker  control  laws  for  the 
KC-135A.  Controllers  are  designed  for  each  of  the  three 
flight  conditions  discussed  in  chapter  III.  This  chapter 
begins  with  a  listing  of  the  requirements  for  the  design 
method,  followed  by  a  detailed  discussion  of  the  maneuvers 
to  be  preformed  and  the  inputs  and  outputs  needed  to  perform 
those  maneuvers.  Next,  the  procedures  used  to  arrive  at  a 
final  design  are  presented  followed  by  the  control  laws 
developed  using  those  procedures.  A  computational  time 
delay  of  one  sampling  period  is  added  and  a  detailed 
comparison  of  the  delayed  results  with  the  undelayed  is 
given.  The  computational  time  delay  results  in  some 
undesirable  effects  such  as  overshoot  in  the  outputs  and/or 
oscillations  in  control  surfaces.  The  controllers  are 
redesigned  so  as  to  minimize  these  undesirable  effects.  In 
order  to  reduce  the  number  of  amplifiers  required  to 
implement  the  designs,  the  controller  matrices  and  K1  are 
studied.  Justification  for  setting  some  of  the  elements  of 
Kq  and  K1  to  zero  is  given  along  with  the  result  of  these 
reduced  controllers.  Finally  robustness  is  tested  by 
applying  the  control  laws  developed  for  flight  condition  #3 
with  computational  time  delay  and  minimum  and  to  each 
of  the  other  two  flight  conditions. 
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Requirements 


The  following  requirements  must  be  satisfied  in 
applying  the  design  method. 

1.  The  number  of  inputs  must  equal  the  number  of 
outputs . 

2.  The  system  must  be  controllable  and  observable. 

3.  All  transmission  zeros  must  be  in  the  left  hand 
s-plane. 


4. 


Rank 


B 

0 


n+p 


Appendix  A  gives  a  state  space  model  for  each  of  the 
three  flight  conditions  studied.  Controllability  and 
observability  of  each  of  the  models  can  be  checked  by  using 
the  controllability  matrix  M  and  the  observability  matrix  M 
respectively.  Where 


Rank  Mc  =  Rank  [B  AB 


n 


(5-1) 


must  be  true  if  the  system  is  controllable  and 

Rank  MQ  =  Rank  [CT  ATCT  ...  AT(n-1)CT]  =  n  (b-2) 


must  be  true  if  the  system  is  observable. 

Transmission  zeros  can  be  checked  using  a  computer 
program  given  by  Reference  7.  Since  the  three  conditions 
studied  all  involved  irregular  designs,  the  matrix  F  = 

[ F2]  must  be  used  in  place  of  C  =  [C 1  C2)  in  checking 
for  the  transmission  zeros. 

Condition  one  above  requires  that  the  number  of  outputs 
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equal  the  number  of  inputs.  This  is  also  true  for  the  three 
flight  conditions  studied  and  a  discussion  of  the  inputs  and 
outputs  is  given  later  in  this  chapter. 

Maneuvers 

Maneuvers  to  be  performed  need  to  meet  two 
requirements.  They  must  be  useful  to  the  pilot  and  they 
must  be  feasible  for  the  aircraft  to  perform.  Remember  that 
this  is  a  heavy  transport  aircraft  not  a  small  highly 
maneuverable  fighter.  In  order  to  illustrate  the 
performance  of  this  design  method,  four  maneuvers,  two 
lateral  and  two  longitudinal,  are  chosen.  The  four 
maneuvers  are: 

1.  Coordinated  turn  (  $  £  30  ,0-0) 

2.  Sideslip  (  |3  <  5,  <t>  =  0) 

3.  Constant  airspeed  climb  (7>  0,  u  =  0) 

4.  Pitch  pointing  (  0  <  5 ,  h  =  0 ,  u  =  0 ) 

The  first  two  maneuvers  are  intended  to  demonstrate  the 
performance  of  lateral  controllers  while  the  next  two 
demonstrate  longitudinal  controller  performance. 

The  coordinated  turn  and  constant  airspeed  climb 
(normal  climb)  are  basic  flight  maneuvers  that  are  performed 
in  every  operational  mission  of  the  KC=135A.  Since  it  is 
desirable  to  have  a  minimum  of  four  maneuvers  to  demonstrate 
the  design  method,  two  more  maneuvers  listed  above  are 
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When  flying  an  approach  to  landing,  under  crosswind 
conditions,  the  KC-135A  aircraft  flies  with  a  crab  into  the 
wind.  Due  to  side  forces  on  the  landing  gear  the  aircraft 
cannot  land  in  this  crab.  Thus  the  pilot  is  required  to  use 
a  wing  low  technique  to  eliminate  the  crab  just  prior  to 
landing.  To  perform  this  the  pilot  must  use  an  appropriate 
amount  of  rudder  to  align  the  aircraft  with  the  runway,  and 
at  the  same  time  bank  the  aircraft  enough  to  prevent  a  side 
velocity  from  developing.  This  maneuver  would  be  very 
useful,  but  it  cannot  be  simulated  with  MULTI  since  it 
requires  nonlinear  sines  of  angles  in  the  output  matrix  C. 
Until  nonlinear  terms  are  allowed  in  MULTI  a  compromise  had 
to  be  made.  Thus  the  choice  of  the  sideslip  maneuver. 
Although  the  sideslip  would  have  limited  use  in  the 
aircraft,  it  does  allow  for  demonstration  of  the  design 
method.  Finally,  during  a  landing  in  which  a  minimum  flap 
setting  is  used,  the  KC-135A  aircraft  attitude  is  such  that 
there  is  the  possibility  of  dragging  the  tail.  Thus  the 
pitch  pointing  is  chosen  in  order  to  alleviate  this 
occurrence . 

Input/Output  Models 

The  "inputs"  of  the  system  refers  to  the  control 
surfaces  of  the  aircraft.  Control  surfaces  for  the  KC-135A 
include  two  lateral  controls,  the  rudder  (6r  )  and  the 

ailerons  and  spoilers  (5^.)  along  with  three  longitudinal 
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controls,  the  elevator  (5  i  ,  the  thrust  (6_),  and  speed 
brakes  ( 6 ) .  These  control  surfaces  are  discussed  in 
detail  in  Chapter  III,  The  designs  developed  in  this 
chapter  include  the  control  surface  limitations  discussed  in 
Chapter  III  and  listed  below. 


-17 

deg  £ 

*r  < 

17 

deg 

(5-3) 

-90 

deg  £ 

'w  i 

90 

deg 

(5-4) 

-25 

deg  £ 

< 

25 

deg 

(5-5) 

0 

deg  £ 

60 

deg 

(5-6) 

4.j.  is  dependent  upon  thrust  available  and  varies  for 
each  flight  condition.  In  general,  the  thrust  required  for 
each  of  the  three  trimmed  flight  conditions  does  not  exceed 
50%  of  available  thrust.  Thus  a  subjective  limit  of  +  50% 
is  placed  on  4^. 

The  "outputs"  are  the  variables  contained  in  the  vector 
y  and  as  such  represent  the  responses  which  are  to  be 
controlled.  The  outputs  can  be  states  of  the  system  or 
linear  combinations  of  the  states.  For  this  study,  the 
lateral  directional  outputs  are  chosen  to  be  the  sideslip 
angle  (0)  and  the  roll  angle  (0) .  Due  to  distinct 
differences  in  the  two  longitudinal  maneuvers  being 
performed,  two  sets  of  longitudinal  outputs  are  chosen.  The 
first,  used  to  perform  the  constant  airspeed  climb,  is  the 
flight  path  angle  (7),  where  y  =  B-af  and  (u),  the 
perturbation  velocity  in  the  x-direction.  The  second  set, 
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used  to  perform  pitch  pointing,  is  the  pitch  angle  (0),  the 
perturbation  velocity  in  the  x-direction,  (u) ,  and  the 
altitude  (h)  where,  for  small  angle  approximations, 
h  =  U  ( d-a ).  Thus  for  the  design  of  the  lateral  controllers 
the  state,  input,  and  output  vectors  are: 


(5-7) 


For  the  longitudinal  case  of  the  constant  airspeed  climb: 


and  for  pitch  pointing: 


(5-8) 


u  =  I  $ 


(5-9) 


Where  e  represents  the  body  bending  at  the  crew  station  as 
discussed  in  Appendix  A. 


Design  Procedures 

What  follows  is  a  brief  discussion  of  the  steps  used  in 
finding  digital  controllers  which  achieve  "good"  performance 


characteristics.  Before  describing  this  it  is  necessary  to 
first  define  the  guide  lines  used  by  this  author  in 

determining  what  "good"  performance  is.  All  command  inputs 
greater  than  zero  are  ramped  up  to  their  steady-state 
values.  This  usually  rosulcs  in  some  overshoot  of  the 
output  before  it  reaches  steady-state.  In  determining 
"good"  responses  it  is  first  desired  to  keep  as  many  outputs 
as  possible  from  overshooting  their  steady-state  values  by 
more  that  2%.  It  is  al.o  desired  to  keep  control  surface 
deflections  to  a  minimum  and  to  prevent  any  unnecessary 
oscillations  in  the  control  surfaces. 

The  steps  used  by  this  author  in  development  of  the 
digital  controllers  presented  in  this  chapter  are  listed 
below.  A  discussion  of  each  step  follows: 

1.  For  each  of  the  state  space  models  given  in 
Appendix  A,  determine  the  design  method  of  Chapter 
II  that  is  to  be  used. 

2.  If  the  design  method  selected  is  irregular,  which 
requires  a  'M'  matrix,  determine  the  minimum  number 
of  M  matrix  elements  required. 

3.  Check  controllability  and  observability  of  the 
system. 

4.  Determine  the  location  of  the  transmission  zeros. 

5.  Select  a  sampling  time  for  discrete  controllers. 

6.  Vary  a,  a,  and  elements  of  M  and  2  to  achieve  a 
"good"  design. 

The  first  step  in  the  design  is  to  determine  which  of 
the  design  methods  discussed  in  Chapter  II  is  to  be  used. 
Since  the  matrices  A,  B,  C,  and  D  of  Equations  (2-1)  and 
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(2-2)  are  known,  the  choice  is  narrowed  to  the  Known/Regular 
or  Known/Irregular  methods.  However,  the  first  Markov 
parameter,  CB,  does  not  have  full  rank,  thus  the 
Known/Irregular  design  method  must  be  used. 

Since  the  design  method  requires  a  measurement  matrix 
'M' ,  the  next  step  is  to  determine  which  elements  of  the  'M* 
matrix  are  to  be  used.  Three  measurement  matrices  are 
selected,  one  for  the  lateral-direction  case,  and  one  each 
for  the  two  maneuvers  in  the  longitudinal  case.  Each  of 
these  'M.'  matrices  are  chosen  based  on  the  requirements 
outlined  in  Chapter  II  and  given  in  Reference  8. 

Next,  controllability  and  observability  are  checked 
using  the  matrices  Mc  and  MQ  as  discussed  earlier.  Also  the 
transmission  zeros  are  determined  using  a  computer  program 
given  by  Reference  7. 

For  the  discrete  systems  studied,  a  sampling  time,  T, 
must  also  be  selected.  The  sampling  time  of  0.05  seconds 
(20  hertz)  is  chosen  since  this  represents  a  reasonable 
sampling  rate  for  aircraft  computers. 

The  final  step  is  to  vary  the  design  parameters  in 
order  to  achieve  a  "good"  design.  Initial  design  are 
accomplished  using  no  actuators  or  sensors.  Using  Options 
11,  13,  12,  and  16  of  the  MULTI  program  the  following  values 
are  entered:  <*=1,  f=l,  all  diagonal  elements  of  2-1,  and  al1 
necessary  elements  of  M=0.25.  Next,  the  matrices  and  K  ^ 
are  examined  for  magnitude.  Generally,  in  order  to  achieve 
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stable  responses  the  relative  magnitudes  of  each  element  in 
these  matrices  should  be  on  the  order  of  10  or  100.  If  the 
magnitudes  are  found  to  be  too  large,  the  value  of  epsilon, 
under  Option  #13,  needs  to  be  reduced.  After  achieving 
stability,  individual  parameters  are  varied  in  order  to 
observe  its  effect  on  system  performance.  The  following 
generalizations  are  found  to  hold  true  under  most 
conditions.  Epsilon,  the  sigma  matrix  multiplier,  affects 
output  overshoot.  Reducing  epsilon  reduces  overshoot. 
Reductions  in  epsilon  usually  require  an  increase  in  alpha, 
which  controls  the  proportion  of  integral  feedback,  in  order 
to  keep  rise  time  of  the  outputs  from  increasing.  Changes 
of  individual  sigma  elements  affect  certain  inputs  and 
outputs.  Increasing  sigma  would  in  general  result  in  an 
increase  in  its  respective  control  surface  deflections  and  a 
reduction  in  rise  time  of  the  outputs.  Increasing  some  M 
matrix  elements  greatly  reduced  control  surface  deflections, 
but  also  slowed  responses.  In  all  design  cases,  notes  are 
taken  on  the  effects  of  vary  each  of  the  design  parameters. 
These  notes  proved  helpful  when  actuators  are  added,  time 
delays  are  considered,  and  when  designing  for  other 
maneuvers  and/or  other  flight  conditions. 

After  a  design  is  found  to  work  with  no  actuators,  the 
next  step  is  to  include  them.  Using  guidelines  found  in 
Reference  12  all  control  surface  actuators  are  entered  as 
10/(s+10)  using  Option  #4.  In  order  to  simulate  engine 
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dynamics  a  transfer  function  of  2/(s+2)  is  added  for  the 
thrust.  In  all  cases,  controllers  developed  with  no 
actuators  produced  unstable  responses  when  actuators  are 
added.  Thus,  it  became  necessary  to  again  trim  the  design 
parameters. 

Next,  the  addition  of  a  computational  time  delay  of  one 
sample  period  is  considered.  In  general,  the  addition  of 
this  time  delay  results  in  two  noticeable  differences. 

These  differences  are  discussed  later  with  examples  being 
presented.  The  final  design  is  then  developed  reducing  any 
ill  effects  caused  by  introduction  of  the  time  delay. 

Tracker  Control  Laws  For  Flight  Condition  #3 

Using  the  above  design  procedure  controllers  are 
developed  for  each  of  the  three  flight  conditions  being 
studied.  A  detailed  discussion  of  the  results  obtained  for 
flight  condition  #3  follows.  This  includes  tabular  listings 
of  Peak  Value,  Final  Value,  tg  and  tp  for  the  outputs  being 
controlled  along  with  Peak  Values  for  each  of  the  control 
surfaces  used.  Also  included  are  Calcomp  plots  showing 
simulation  results. 

Coordinated  Turn ,  Flight  Condition  #3  To  accomplish  a 
coordinated  turn,  the  bank  angle,  phi,  is  commanded  to  30 
degrees  in  4  seconds  while  sideslip,  beta,  is  commanded  to 
zero.  The  control  law  selected  for  this  maneuver  is  defined 
by  the  following  parameters: 
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T  = 
a  = 

t  = 
2  = 

M  = 


0.05  seconds 

(20  hertz) 

(5-10) 

4.25 

(5-11) 

0.045 

(5-12) 

diagonal  [1 

.  5 

2.0] 

(5-13) 

0.75  0.0 

(5-14) 

0.0  0.4 

which  yield 


0.02204 
0.9891 
.005186 
.2327 


1.68 

-1.486 

0.3953 

-0.3496 


(5-15) 

(5-16) 


Figures  5-la  thru  5-lf  show  the  simulation  for  this 
controller  for  a  commanded  30  degree  coordinated  turn. 

The  command  input  vector  is  given  by  Figure  5-la. 

Figure  5-2b  shows  the  roll  angle,  phi.  The  peak  value 
30.478  degrees  occurs  at  7  seconds.  This  is  3  seconds  after 
the  commanded  input  reaches  steady-state.  Settling  time  for 
phi  is  5.2  seconds  and  is  based  on  2%  of  the  final  value. 
Sideslip,  shown  in  Figure  5-lc,  is  almost  negligible  with 
peak  values  of  +0.015-  and  -0.04.  Figures  5-le  and  5-ld  show 
control  wheel  and  rudder  deflections  respectively.  The 
control  wheel  proved  to  be  a  big  factor  in  performing  this 
maneuver.  Settling  time  for  phi  had  to  be  compromised  in 
order  to  keep  maximum  deflections  of  the  control  wheel 
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within  acceptable  limits.  A  reminder  of  what  the  control 
wheel  represents  is  needed  at  this  time  for  clarification  of 
what  Figure  5-ld  represents.  The  peak  value  of 
approximately  30  degrees,  in  Figure  5-ld,  represents  1/3  of 
total  control  wheel  deflection.  Assuming  a  linear 
representation  between  the  control  wheel  and  aileron 
movement,  this  would  represent  approximately  7  degrees  of 
aileron  deflection.  This  7  degrees  of  movement  occurs  in  0.3 
seconds,  which  equates  to  a  rate  of  23.3  deg/sec-well  within 
the  capability  of  the  actuator.  Likewise,  a  rudder 
deflection  of  6  degrees  in  3.8  seconds,  as  shown  in  Figure 
5-le,  equates  to  a  deflection  rate  of  1.57  deg/sec.  Again, 
well  within  the  limits  of  the  actuator. 

Additionally,  based  on  the  assumptions  used  in  Appendix 
A,  Blakelock  (Ref.  12)  requires  that  a  constant  yaw  rate,  r, 
be  achieved  in  a  coordinated  turn.  This  yaw  rate  is  defined 
by: 

r  =  57 .3  [g/Ug]  sin  <P  deg/sec  (5-17) 

Figure  5-lf  shows  that  the  state  r  does,  in  fact,  reach 
the  steady-state  value  of  3.925  deg/sec  as  computed  from 
Equation  (5-17).  Table  5-1  summarizes  the  results  of  the 
simulation  represented  by  Figures  5-la  thru  5-lf. 
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TABLE  5-1 

Simulation  Results  For 
Coordinated  Turn  Flight  Condition  #3. 


Input/ 

Output 

Peak 

Value 

F  inal 
Value 

tp( sec) 

t  (sec) 

<t>  (deg) 

30.478 

30.0 

7 

5.2 

0  (deg) 

-0.04 

-0.032 

7 

14.4 

r  (deg/sec) 

3.93 

3.93 

5.2 

5.2 

Sr  (deg) 

-5.5 

— 

— 

— 

ydeg) 

32.0 

-- 

— 

-- 

Sideslip,  Flight  Condition  #3 .  The  control  laws 
defined  bv  Equations  (5-10)  thru  (5-16)  are  used  to  simulate 
the  sideslip  maneuver.  This  maneuver  is  accomplished  by 
commanding  sideslip  to  5  degrees  in  8  seconds  while 
commanding  bank  angle  to  zero.  When  this  maneuver  is 
performed,  the  bank  angle  exceeds  two  degrees.  In  order  to 
improve  this  result  a  new  controller  is  designed,  a 
constraint  imposed  on  this  new  controller  is  to  keep  the 
number  of  gain  changes  in  H,  K  ,  and  K  given  by  Equations 
( 5-14 ),( 5-15 ) ,  and  (5-16)  to  a  minimum.  Thus,  the  following 
parameters  are  used  in  design  of  the  new  controller: 


T  = 

0.05  seconds  (20  hertz) 

(5-18) 

a  = 

4.25 

(5-19) 

€  = 

0.045 

(5-20) 

v  _ 

<U 

d iag [1.5  2.0] 

(5-21) 

►  *  *  1 
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0.25  0.0 


(5-22) 


(5-23) 

(5-24) 


Notice  that  this  new  controller  involves  only  changes 
in  the  first  columns  of  M,  Kq,  and  Kj .  Using  this 
controller  design  the  sideslip  maneuver  is  accomplished  with 
the  results  of  the  simulation  shown  in  Figures  5-2a  thru 
5-2d . 

The  outputs  for  this  maneuver  are  shown  in  Figure  5-2b 
for  the  inputs  of  Figure  5-2a.  The  sideslip  angle  peaked  at 
5.076  degrees  at  10  seconds  and  settles  to  within  2%  of  the 
final  value  of  5  degrees  by  8.4  seconds,  0.4  seconds  after 
the  input  reaches  steady-state.  Roll  angle  is  negligible 
with  a  peak  of  only  -0.2426  occurring  at  8.2  seconds. 

The  control  surfaces  are  shown  in  Figure  5-2c.  Both 
have  initial  transients  and  transients  at  8  seconds 
corresponding  to  slope  changes  in  the  command  input.  These 
transients  are  considered  minimal.  The  transients  occurring 
at  8  seconds  can  be  reduced,  but  at  the  expense  of 
increasing  settling  time.  Therefore  a  compromise  must  be 
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FIGURE  5-2a : 


TIME.  SECONDS 


Flight  Condition  #3,  Sideslip 
(Command  Input  Vector) 


I  ME.  SECONOS 


FIGURE  5-2 b : 


Flight  Condition  #3,  Sideslip 
(Output  1:  Roll  Angle;  Output  2;  Sideslip) 
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made  here 
Table  5-2, 


The  resu.  of  this  simulation  are  summarized  in 


Table  5-2 

Simulation  Results  For 
Sideslip  Flight  Condition  #3. 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec 
P 

<P  (deg) 

-0.242 

0.0 

8.2 

0  (deg) 

5.076 

5.0 

10.0 

K l deg ) 

7.8 

— 

— 

6w(deg) 

28.0 

— 

— 

t  (sec) 
s 


Normal  Climb,  Flight  Condition  JCK  Using  the  design 
procedures  discussed  earlier,  a  longitudinal  controller  is 
designed  to  accomplish  a  1200  fpm  climb.  To  accomplish  this 
flight  'ath  angle,  gamma,  is  commanded  to  a  value  of  4.86 
degrees  in  b  seconds  and  at  the  same  time  perturbation 
velocity  in  the  x-direction  is  commanded  to  zero.  The 
control  law  selected  for  this  maneuver  is  defined  by  the 
parameters: 


0.05  seconds  (20  hertz) 
2.5 
0 , 5 


> 

2  =  d  lag [0.5 

0.3] 

b\ 

- 1 

O 

• 

O 

G.75 

O 

• 

O 

r- 

a 

K  =1 

f \\ 

l 

O 

« 

O 

O 

O 

O 

* 

O 

(5-25) 

(5-26) 

(5-27) 

(5-28) 

(5-29) 


as 


which  yields 


-0.9558  0.0 

0.04861  2.957 


’-0.3823  0.0 

0.01944  1.183 


(5-30) 

(5-31) 


Figures  5-3a  thru  5-3h  show  the  simulation  results  for 
this  controller.  The  command  input  vector  is  given  by 
Figure  5-3a.  Figure  5-3b  shows  flight  path  angle  and 
velocity.  Flight  path  angle  has  a  peak  value  of  5.32 
occurring  at  8.1  seconds  which  is  2.1  seconds  after  the 
input  reaches  steady-state.  Settling  time  for  the  flight 
path  angle  is  12.3  seconds/  6.3  seconds  after  the  input 
reaches  steady-state.  At  first  this  may  appear  to  be  a 
rather  large  settling  time.  But  when  the  size  of  the 
aircraft  and  the  slow  responding  engines  are  considered, 
12.3  seconds  to  stabilize  in  a  1200  rpm  constant  airspeed 
climb  isn't  unreasonable.  Perturbation  velocity,  u,  peaked 
at  -0.1476  ft/sec  in  4.8  seconds,  then  slowly  returns  to 
zero. 

The  flight  controls  needed  for  this  maneuver  are  shown 
in  Figure  5-3c.  Thrust  peaks  at  24%  in  7  seconds.  This 
should  be  easily  attainable  by  the  engines  on  the  KC-135A. 
Remember  that  this  represents  an  increase  in  thrust  above 
that  used  for  level  flight.  The  elevator  averages  a  3 
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FIGURE  5-3b:  Flight  Condition  #3,  Normal  Climb 

(Output  1:  Flight  Path  Angle;  Output  2:  Veloc 


degree  deflection  over  a  5  second  duration  and  the  rates  of 
deflection  are  approximately  20  deg/sec.  Again,  an  easily 
attainable  limit. 

Figure  5-3d  shows  the  altitude  h  as  a  function  of  time. 
The  slope  of  this  curve  represents  the  desired  1200  fpm 
climb.  Figures  5-3f  thru  5-3h  represent:  the  body  bending 
mode.  This  body  bending  is  measured  at  the  crew  station. 
Figure  5-3f  shows  the  first  bending  mode,  Figure  5-3g  shows 
the  second  and  Figure  5-3h  shows  the  sum  of  the  two. 
Comparison  of  Figure  5  —  3  f  and  Figure  5-3h  shows  that  the  1st 
bending  mode  accurately  represents  the  total  bending  of  the 
aircraft  and  that  future  studies  involving  aircraft  bending 
need  only  consider  this  mode.  This  by  no  means  represents 
any  new  discovery  but  merely  supports  what  has  already  been 
known.  It  should  also  be  noted  that  body  bending  is  only 
observed  here  and  that  its  effects  on  sensors,  etc.  are  not 
included  in  the  model.  See  Appendix  A  and  Reference  10  for 
the  equations  used  to  model  the  aircraft  bending.  Table  5-3 
gives  a  summary  of  the  results  of  this  simulation. 


Table  5-3 

Simulation  Results  For 
Normal  Climb  Flight  Condition  #3. 
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Pitch  Pointing ,  Flight  Condition  #3 .  The  longitudinal 
controller  used  to  perform  the  normal  climb  uses  only 
elevator  and  thrust  as  inputs.  Speed  Brakes  are  left  out  of 
the  model  since  they  would  contribute  nothing  to  this 
maneuver.  For  pitch  pointing,  speed  brakes  need  to  be 
included  since  there  are  now  three  outputs  to  control  and 
thus  three  inputs  are  needed.  Therefore,  a  second 
longitudinal  controller  is  designed  using  three  inputs. 

This  controller  is  defined  by  the  parameters: 


T  =  0.05  seconds  (20  hertz)  (5-32) 

«  -  2.5  (5-33) 

C  =  0.3  (5-34) 

2=  d iag [1.0  0.75  0.1]  (5-35) 

0.75  0.0 

M  =  0.0  0.25  (5-36) 

0.0  0.0 
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vcv.v. 


4.V.V^ 


-0.4157  -2.597  0.0 

K  =  -4.017  -1.508  0.0  (5-37) 

u  -1.978  -0.6188  0.5915 

» 

-0.1663  -1.039  0.0 

K.  =  -1.607  -0.6033  0.0  (5-38) 

-0.7914  -0.2475  0.2366 

- 

To  perform  pitch  pointing,  pitch  angle  is  commanded  to 
4  degrees  in  4  seconds  while  at  the  same  time  altitude  and 
perturbation  velocity  are  commanded  to  zero.  Figures  5-4a 
thru  5-4h  represent  the  simulation  results  for  this  maneuver 
using  the  controller  defined  by  Equations  (5-32)  thru 
(5-38).  The  command  input  vector  is  given  in  Figure  5-4a. 
Figure  5-4b  shows  two  of  the  outputs,  altitude  and 
perturbation  velocity.  Both  of  these  remain  negligible 
since  altitude  peaks  at  0.299  feet  in  4.4  seconds  and  u 
peaks  at  -0.479  in  1.6  seconds.  Figure  5-^0  shows  pitch  and 
attack  angles.  Pitch  angle  closely  follows  the  input  and 
reaches  a  peak  of  4  degrees  in  16  seconds.  Settling  time 
for  pitch  angle  is  4.4  seconds.  Elevator  and  thrust  are 
shown  in  Figure  5-4d.  Thrust  peaks  at  32  %rpm  in 
approximately  4.5  seconds.  The  elevator  is  somewhat  shaky 
and  thus  decreases  the  effectiveness  of  this  maneuver. 

Figure  5-4e  shows  speed  brake  deflections  which  peak  at  52 
degrees  in  4.25  seconds.  This  in  itself  says  something 
about  the  maneuver  since  53  degrees  of  speed  brakes  to 
produce  4  degrees  of  pitch  pointing  doesn't  seem  like  much 
of  a  trade  off.  Figures  5-4f  thru  5-4h  show  the  body 
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Table  5-4  summarizes 


VN 


bending  modeled  at  the  crew  station, 
the  results  of  this  maneuver. 


Table  5-4 

Simulation  Results  For 
Pitch  Pointing  Flight  Condition 


#3 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec) 

P 

t  (sec) 
s 

h  (ft) 

0.299 

0.0 

4.4 

— 

6  (deg) 

4.0 

4.0 

16.0 

u  (ft/sec) 

-0.479 

0 . 0 

1.6 

— 

5e  (deg) 

4.0 

— 

— 

— 

Wdes) 

52.0 

— 

— 

— 

5t(%RPM) 

32.0 

-- 

-- 

— 

Tracker  Control  Laws  For  Flight  Condition  #2  and  #1 

Using  the  same  procedures,  controllers  are  developed 
for  each  of  the  two  remaining  flight  conditions.  What 
follows  is  a  list  of  the  parameters  used  to  define  each  of 
the  controllers  used.  Simulations  are  accomplished  for  each 
of  these  controllers.  Rather  than  discuss  these  simulations 
is  detail,  their  results  are  simply  compiled  in  Appendix  C 
for  observation. 

Coord ina ted  Turn ,  Flight  Condition  #2  The  control  law 
developed  for  this  maneuver  is  defined  by  the  following 
parameters : 


V3 
*  - 

VV-' 
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•  V 


T 


0.05  seconds  (20  hertz) 


(5-39) 


’ 
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Figures  C-2a  thru  C-2d  and  Table  C-2  given  in  Appendix 
C  show  the  simulation  results  for  this  controller. 

Normal  Climb,  Flight  Condition  #2  The  control  law 
developed  for  this  maneuver  is  defined  by  the  following 
parameters: 


T  =  0.05  seconds  (20  hertz) 
a  =  2.5 
t  =  0.3 

I  =  diagonal  [1.75  .1] 


0.0  0.75 

O 

• 

o 

M  = 

0.0  0.0 

0.0 

o 

n 

-1.01 

o 

• 

o 

0.3118 

2.206 

-0.4039 

o 

• 

o 

K1  = 

0.1247 

0.8824 

Figures 

:-3a  thru 

C-3h  a 

(5-53) 

(5-54) 

(5-55) 

(5-56) 

(5-57) 

(5-58) 

(5-59) 


C  show  the  simulation  results  for  this  controller. 

P itch  Pointing ,  Flight  Condition  #2  The  control  law 
designed  for  this  maneuver  is  defined  by  the  following 
parameters : 


0.05  seconds  (20  hertz) 

(5-60) 

2.5 

(5-61) 

0.3 

(5-62) 

vr-.“  v. 


:  -r.  *\  wi  w  w  •*-.  tx  T*  ~i  “"v  '"X.tv '•■^Ti'nv.T-_''0'V .  ;  *v*r  v- 


irftlWMftBi 


£  =  diagonal  [1  0.75  0.1] 


0.75 

0.0 

0.0 


.0 

.25 


■0.3687 

■2.973 

■6.865 

•0.1475 

■1.189 

■2.746 


-1.398 
-G  .9532 
-1.806 


0.0 

0.0 

2.206 


-0.5592  0.0 

-0.3813  0.0 

-0.7222  0.8824 


(5-63) 


(5-64) 


(5-65) 


(5-66) 


Figures  C-4a  thru  C-4h  and  Table  C-4  given  in  Appendix 
C  show  tne  simulation  results  for  this  controller. 

Coordinated  Turn ,  Flight  Condition  #1  The  control  law 
developed  for  this  maneuver  is  defined  by  the  following 
parameters : 


=  0.05  seconds  (20  hertz) 


=  0.03 


=  diagonal  [1.5  5.0] 


0.75  0.0 

0.0  0.4 


0.00525  1.488 


0.6493  -2.527 


0.001312  0.3721 


0.1623  -0.6317 


(5-67) 

(5-68) 

(5-69) 

(5-70) 


(5-71) 


(5-72) 


(5-73) 


Figures  C-5a  thru  C-5f  and  Table  C-5  given  in  Appendix 
C  show  the  simulation  results  for  this  controller. 


'■.•'■./O'jOvOv 
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S  idesl  ip,  Flight  Cond  i  t  ion  £1  The  control  law 
developed  for  this  maneuver  is  defined  by  the  following 


parameters : 

T  -  0.05  seconds  (20  hertz)  (5-74) 

a  =  4.0  (5-75) 

t  =  0.03  (5-76) 

2  =  diagonal  [1.5  5]  (5-77) 

0.25  0.0 

M  =  (5-78) 

0.0  0.4 


0.01575 

1.948 


1.488 
-2 . 527 


“0.003937  0.3721 

0.487  -0.6317 


(5-79) 

(5-80) 


Figures  C-6a  thru  C-6d  and  Table  C-6  given  in  Appendix 
C  show  the  simulation  results  for  this  controller. 

Normal  Climb,  Flight  Condition  #1  The  control  law 
developed  for  this  maneuver  is  defined  by  the  following 
parameters: 


T  =  0.05  seconds  (20  hertz) 
a  =  2.5 


t  =  0.3 

2  =  diagonal  [1.5  0.2] 


M 


0.0  0.75  0.0 

0.0  0.0  0.0 


-0.9256 

0.2171 


0.0 

3.028 


(5-81) 

(5-82) 

(5-83) 

(5-84) 

(5-85) 

(5-86) 
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1-0.3702  0.0 

K,  =  (5-87) 

I  C .08682  1.211 

Figures  C-7a  thru  C-7h  and  Table  C-7  given  in  Appendix 
C  show  the  simulation  results  for  this  controller. 

P itch  Pointing ,  Flight  Condition  #1  The  control  law 
developed  for  this  maneuver  is  defined  by  the  following 
parameters : 


T  =  0.05  seconds  (20  hertz)  (5-88) 

<x  s  4.5  (5-89) 

<  =  0.3  (5-90) 

-  diagonal  [1  0.75  0.1]  (5-91) 

0.75  0.0 

M  =  0.0  0.25  (5-92) 

0.0  0  o  0 

-0.272  -1.467  0.0 

Kn  =  -2.588  -0.9788  0.0  (5-93) 

-4.257  -1.29  1.514 

-0.1088  -0.5869  0.0 

K.  =  -1.035  -0.3915  0.0  (5-94) 

-1.703  -0.516  0.6056 


Figures  C-8a  thru  C-8h  and  Table  C-8  given  in  Appendix 
C  show  the  simulation  results  for  this  controller. 


Time  Dela\ 


Once  the  control  laws  are  developed  for  each  of  the 
maneuvers  the  next  step  in  achieving  a  final  design  is  to 
consider  a  computational  time  delay  of  one  sampling  period. 
Computational  time  delay  represents  the  aircraft  computer 
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computation  time  needed  to  generate  the  input  to  the 
actuators,  based  on  the  current  error  vector.  No  delay 
would  mean  that  the  computer  performs  instantaneous 
calculations  of  the  control  law.  A  delay  of  one  sampling 
period  means  that  the  time  it  takes  the  computer  to  perform 
these  calculations  is  less  than  one  sampling  time.  This 
time  delay  is  entered  under  Option  #23  in  MULTI.  What 
follows  is  a  discussion  of  the  effect  of  uhe  time  delay  for 
each  of  the  previously  designed  controllers  of  flight 
condition  #3. 

C-oordina -ed  Turn,  Flight  Condition  #3  With  Time  Delay ♦ 
Using  the  control  law  given  by  Equations  (5-10)  thru  (5-16) 
the  coordinated  turn  is  performed  again,  this  time  including 
a  computational  time  delay  of  one  sampling  period.  Figures 
5-5a  thru  5-5f  shows  the  simulation  results  for  this 
condition.  Comparison  of  Figures  5-5a  thru  5-5f  with 
Figures  5-la  thru  5-lf  show  the  effects  of  the  computational 
time  delay.  At  first  glance  there  appears  to  be  very  little 
difference,  but  a  comparison  of  the  figures  of  merit  and  the 
maximum  control  deflections  shows  that  there  are  indeed  some 
differences.  Tabje  5-5  gives  a  summary  of  the  results  of 
the  computational  time  delay  simulation. 
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Table  5-5 

Simulation  Results  For 

Coordinated  Turn  With  Time  Delay  Fliyht  Condition  #3. 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec) 

P 

t  (sec) 

<t>  (deg) 

31.229 

30.0 

6.6 

10.0 

0  (deg) 

0.099 

-0.03 

2.2 

15.2 

r  (deg/sec) 

4.10 

3.925 

5.2 

10.0 

4r(deg) 

-5.5 

— 

— 

*w(deg) 

32.0 

— 

— 

— 

By  comparison  of  Table  5-5  with  Table  5-1  the  effects 
of  the  time  delay  become  more  noticeable.  from  this  it  is 
seen  that  the  bank  angle  has  an  overshoot  of  1.229  degrees 
and  that  settling  time  has  increased  to  10  seconds.  Control 
deflections  remain  about  the  same  and  sideslip  is  still 
negligible.  These  effects  appear  to  be  acceptable  and 
nothing  mere  needs  to  be  done,  but  if  they  are  not,  then,  it 
is  necessary  to  redesign  the  controller  to  reduce  or 
eliminate  any  of  the  undesired  time  delay  effects.  Either 
way,  the  overshoot  of  the  bank  angle  can  be  controlled  with 
the  development  of  a  new  controller.  To  demonstrate  that 
this  is  possible,  using  the  same  guidelines  as  before  a  new 
controller  is  designed  using  the  following  parameters; 


T  =*  0.05  seconds  (20hertz) 
a  =  4.25 


(5-95) 

(5-96) 


(  =  0.08 


M 


K, 


diagonal  [1.5  2.0] 


0.75  0.0 

0.0  1.0 

*0.03104 

1.337 

1.766 

-1.182 

0.007304 

0.3145 

0.4154 

-0 . 2782 

(5-97) 

(5-98) 

(5-99) 

(5-100 ) 

(5-101) 


Figures  5-0  thru  5-6f  show  the  simulation  results  for 
this  new  controller.  The  biggest  problem  encountered  in 
trying  to  reduce  the  bank  angle  overshoot  caused  by  the 
computational  time  delay  is  keeping  control  surfaces  from 
oscillating.  Table  5-6  summarizes  thi  results  of  the  rede¬ 
signed  controller  given  by  Equations  (5-95)  thru  (5-101). 


Table  5-6 

Simulation  Results  For  Coordinated  Turn  With 
Time  Delay  Flight  Condit'on  #J  (Redesigned) 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec) 

P 

t  (sec) 
s 

<#>  (deg) 

30 . 563 

30.0 

7.0 

5.2 

0  (deg) 

O 

« 

o 

to 

-0.082 

8.4 

14.0 

r  (deg/sec) 

4.17 

3.93 

5 . 0 

7.2 

(deg) 

-5.5 

— 

— 

— 

5w(deg) 

32.0 

— 

— 

(0EGRLLS1 


FIGURE  5~6a : 


Flight  Condition  #3,  Coordinated  Turn 
With  Time  Delay,  Redesigned 
(Command  Input  Vector) 


GREES/SEC }  m  CONTROL  WHEEL  ( OEOREES ) 


A  Comparison  of  Table  5-1  with  5-5,  and  5-6  show  that 
with  the  new  controller  the  time  delayed  system's 
performance  compares  vary  favorably  with  the  undeiayed 
system,  thus  showing  that  any  ill  effects  caused  by  the 
introduction  of  the  computational  time  delay  can  be 
eliminated . 

S ideslip.  Flight  Condition  #3  With  T ime  Delay  Using 
the  control  law  given  by  Equations  (5-18)  thru  (5-24)  the 
sideslip  maneuver  is  performed  with  a  controller 
computational  time  delay  of  one  sampling  period.  Figures 
5-7a  thru  5-7d  show  the  results  of  this  simulation.  Also, 
Table  5-7  summarizes  these  results. 
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Table  5-7 

Simulation  Results  For 

Sideslip  With  Time  Delay  Flight  Condition  #3 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  ( sec 
p 

<f>  (deg) 

-0.4879 

0.0 

8.4 

0  (deg) 

5.0359 

5.0 

16.0 

\  (deg) 

8.0 

— 

— 

Sw(deg) 

30.0 

— 

— 

ts (sec) 


.In  this  case  a  comparison  of  Tables  5-2  with  5-7  shows 
that  performance  is  not  lost  when  the  time  delay  is 
included..  However,  by  looking  at  Figure  5-2c  and  5-7c  it  is 
seen  that  the  control  surfaces  tend  to  be  more  oscillatory 
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FIGURE  5-7c: 


T [ME.  SECONDS 


Flight  Condition  #3,  Sideslip  With  Time  Delay 
(Control  1:  Rudder;  Control  2:  Control  Wheel) 


FIGURE  5-7d : 


Flight  Condition  #3,  Sideslip  With  Time  Delay 
(State  3:  Roll  Rate?  State  4:  Yaw  Rate) 
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when  the  time  delay  is  introduced.  Thus  a  new  controller  is 
designed  with  these  oscillations  in  mind.  The  new 
controller  is  defined  by  the  parameters: 

T  = 
at  = 

€  = 

2  = 

M  = 

K0  - 

K,  = 


0.05  seconds 

(20  hertz)  (5- 

102) 

3.0 

(5- 

103) 

0.08 

(5- 

104) 

d iagonal 

[1 

1.2] 

(5- 

105) 

0.25  0 

.0 

(5- 

106) 

O 

o 

• 

o 

.  5 

0.05176 

1 

«• 

.049 

(5- 

107) 

2.486 

-0 

.8270 

0.01725 

0 

.3497 

(5- 

•108) 

0.8285 

m 

-0 

.3093 

5-8a  thru 

5- 

-8d  show  the  simulation  results 

of 

this  new  controller.  A  comparison  of  Figure  5-7c  and  5-8c 
shows  that  it  is  possible  to  reduce  the  oscillations  to  some 
degree.  Table  5-8  summarizes  the  results  for  this  new 
controller.  This  table  shows  that  the  reduction  of  the 
control  surface  oscillations  is  at  the  expense  of  a  slight 
increase  in  settling  time.  These  oscillations  may  be 
reduced  more,  however,  settling  time  will  increase. 
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Flight  Condition  #3,  Sideslip  With  Time  Delay 
Redesigned  (Output  1:  Roll  Angle;  Output  2: 

S  idesl ip) 


■V  >  '-'A  N1'.' 
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ROLL  4  row  RATES  ( OEGREES/SEC )  £  RUOOER  t  CONTROL  WHEEL  DECREES 

.20  -0.80  -0.40  0.00  0.40  0.80  2  -.-8.00  O.CO  8. CO  16.00  24.  CC 


Table  5-R 

Simulation  Results  For  Sideslip  with 
Time  Delay  Flight  Condition  #3  (Redesigned) 


Input/ 

Output 

Peak 

Value 

F  inal 
Value 

t  (sec) 

t  (sec) 

O 

<t>  (deg) 

-0.4475 

0.0 

8.4 

— 

0  (deg) 

5.088 

5.0 

16.0 

8.8 

5r  (deg) 

6.0 

— 

— 

— 

5w(deg) 

28.0 

— 

— 

— 

Normal  Climb,  Flioht  Condition  #3  With  Time  Delay 
Using  the  control  law  defined  by  Equations  (5-25)  thru 
(5-31)  the  normal  climb  with  a  computational  time  delay  of 
one  sampling  period  is  simulated.  Figures  5-9a  thru  5-9h 
show  the  results  sf  this  simulation  and  Table  5-9  summarizes 
the  results. 


Table  5-9 

Simulation  Results  For  Normal  Climb 
With  Time  Delay  Flight  Condition  #3 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec) 

P 

t  (sec) 

7  (deg) 

5.434 

4.86 

8.1 

12.9 

u  (ft/sec) 

-0.3108 

0.0 

4.8 

— 

6e  (deg) 

-3 . 5 

— 

— 

-- 

<5t'%RPM) 

24.0 

— 

— 

— 

VECTOR 


:2.co  ip.cc  z\ 
TIME.  SECONDS 


FIGURE  5-9c:  Flight  Condition  #3,  Normal  Climb  With  Time 

Delay,  (Control  1:  Elevator;  Control  2;  Thrust) 


12. CO  16  CO  20. CC 

H  HE .  SECONDS 


FIGURE  5-9d: 


Flight  Condition  #3,  Normal  Climb  With  Time 
Delay,  (State  1:  Altitude) 
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HME,  SECONDS 


FIGURE  5-9g: 


Flight  Condition  #3,  Normal  Climb  With  Time 
Delay,  (State  9:  2nd  Body  Bending  Mode) 


12. CO  16. CO  20. CO 

n  ME .  SECONDS 


FIGURE  5 - 9 h 


Flight  Condition  #3,  Normal  Climb  With  Time 
Delay,  (Sum  of  1st  &  2nd  Bending  Mode) 


Comparison  of  Table  5-9  with  Table  5-3  and  Figures  5-9a 


thru  5-9h  with  Figures  5-3a  thru  5-3h  show  that  only  minor 
changes  occur  when  introducing  the  time  delay.  Since  these 
changes  are  minor  in  nature  there  is  no  need  to  redesign  the 
controller . 

Pitch  Pointing ,  Flight  Cond  it.  ion  #3  With  Time  Delay 
Using  the  control  law  defined  by  Equations  (5-32)  thru 
(5-38)  the  pitch  pointing  maneuver  including  a  computational 
time  delay  is  performed.  Figures  j-10a  thru  5-L0h  and  Table 
5-10  show  the  rest  _s  of  this  simulation. 

Comparison  of  Table  5-10  with  Table  5-4  and  Figures 
5-10a  thru  5-10h  with  Figures  5-4a  thru  5-4h  shows  the 
effect  of  the  computational  time  delay  on  this  maneuver. 

The  only  substantial  difference  is  in  the  settling  time  for 
the  pitch  angle  which  increases  by  0.8  seconds. 


Table  5-10 

Simulation  Results  For 

Pitch  Pointing  With  Time  Delay  Flight  condition  #3 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec) 

t „ (sec) 

h  (ft) 

0.598 

0.0 

4.6 

«p» 

e  (deg) 

4.0 

4.0 

16.0 

5.: 

u  (ft/sec) 

-0.899 

0.0 

3.0 

— 

4e  (deg) 

5.0 

— 

— 

— 

isb(deg) 

52.0 

— ■ 

— 

— 

5T(%RPM) 

32.0 

—  — 

— 

— 
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FIGURE  5-IOd:  Flight  Condition  #3,  Pitch  Pointing 

With  Time  Delay  (Output  1:  Altitude; 
Output  2:  Velocity) 


FIGURE  5-10d 


Flight  Condition  #3,  Pitch  Pointing 
With  Time  Delay  (Control  1:  Elevator; 
Control  2;  Thrust) 
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FIGURE  5-10g : 


Flight  Condition  #3,  Fitch  Pointing 
With  Time  Delay  (State  9:  2nd  Body 
Bending  Mode) 


FIGURE  5-10h: 


SlCCNDS 


i light  Condition  #3,  Pitch  Pointing 
With  Time  Delay  (Sum  of  1st  &  2nd 
Bending  Modes) 


A  new  controller  is  designed  attempting  to  reduce  this 
settling  time.  As  in  the  other  maneuvers,  when  trying  to 
redesign  the  controller  to  reduce  the  effects  of  the  time 
delay  the  control  surfaces  prove  to  be  very  sensitive.  A 
small  change  in  design  parameters  would  in  this  case  cause 
unacceptable  elevator  and  speed  brake  oscillations.  Thus  a 
compromise  between  the  settling  time  and  control  surface 
movement  is  needed.  If  settling  time  is  allowed  to  increase 
slightly,  elevator  and  speed  brake  oscillations  are  reduced. 
A  new  controller  is  designed  to  show  this  effect.  This 
controller  is  defined  by  the  following  parameters: 


results  of  this  new  controller 
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SURE  5-llc : 


Flight  Condition  #3,  Pitch  Pointing 
With  Time  Delay,  Redesigned 

(State  2:  Pitch  Angle;  State  4:  Attack  Angle) 


FIGURE  5-lld 


TIME.  SECONDS 


Flight  Condition  #3,  Pitch  Pointing 
With  Time  Delay  Redesigned 
(Control  1:  Elevator;  Control  2:  Thrust) 
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(Control  3:  Speed  Brakes) 
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FIGURE  5-llg 


Flight  Condition  #3,  Pitch  Pointing 
With  Time  Delay,  Redesigned 
(State  9:  2nd  Body  Bending  Mode) 
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FIGURE  5-llh : 


Flight  Condition  #3,  Pitch  Pointing 
With  Time  Delay,  Redesigned 
(Sum  of  1st  &  2nd  Bending  Modes) 
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Taole  5-11 

Simulation  Results  For  Pitch 
With  Time  Delay  Flight  Condition 


Point ing 
#3  Redesigned 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec) 

ts (sec) 

h  (ft) 

0.343 

0.0 

4.6 

— 

9  (deg) 

4.0 

4.0 

16.0 

5.4 

u  (ft/sec) 

-0.485 

0.0 

2.4 

— 

*e(deg) 

3.0 

— 

— 

— 

4sl(deg) 

48.0 

— 

— 

— 

*t(%rpm) 

31.0 

—  — 

— 

Comparison  of  Table  5-11  with  5-10  and  Figures  5-lla 
thru  5-llh  with  5-10a  thru  5-llh  demonstrates  that  a  small 
increase  in  settling  time,  0,2  seconds  to  be  exact,  can 
reduce  control  surface  deflections.  Further  increase  in  the 
settling  time,  to  6.5  seconds,  reduces  the  elevator 
deflections  shown  in  Figure  5-lld  even  further.  When 
settling  time  is  allowed  to  increase  to  6.5  seconds,  the 
elevator  has  just  two  spikes  and  the  transients  that  appear 
after  the  two  spikes  shown  in  Figure  5-lld  disappear 
completely.  Thus  a  trade  off  between  control  surface 
movement  and  maneuver  performance  is  required.  This  in 
general  is  found  to  hold  true  for  all  the  resigns  performed 


in  this  thesis 


Mod i f icat ions  to  Kg  and  K ^ 

MULTI  has  been  modified  recently  to  allow  changing  of 
any  element  in  any  matrix  at  any  time  without  having  to 
reenter  t  ,e  entire  matrix.  This  was  originally  added  to 
avoid  the  inconvenience  involved  in  making  an  error  when 
entering  the  A,  B,  and  C  matrices  in  Option  #3.  However, 
this  modification  also  added  an  interesting  sidelight.  The 
matrices  and  Kj  are  stored  in  MULTI  and  can  be  viewed  by 
selecting  Option  #114.  Allowing  the  user  the  option  of 
changing  the  elements  of  these  matrices  in  Option  #114,  the 
effect  of  each  element  on  the  simulation  results  can  be 
tested.  In  fact,  a  reduction  in  amplifiers  needed  to 
implement  the  controller  would  result  if  some  of  these 
elements  could  be  set  to  zero,  without  affecting  the  results 
significantly. 

Examination  of  the  theoretical  basis  of  the  elements  of 
Kq  and  Kj  gives  some  interesting  insight  into  selecting 
elements  to  be  set  to  zero.  To  demonstrate  this,  take  the 
case  of  Equations  (5-100)  and  (5-101)  which  are  repeated 
here  for  convenience. 


for 


0.03104  1.337 

1.766  -1.182 

. 

0.007304  0.3145 


(5-116) 


[0.4154  -0.2782J 

Remember  that  these  represen 
the  coordinated  turn  maneuver 


(5-117) 

t  the  Kg  and  Kj  matrices 
which  was  redesigned  to 
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account  for  the  time  delay.  The  first  column  of  KQ 
represents  gains  applied  to  the  rudder,  while  the  second 
column  represents  gains  applied  to  the  control  wheel.  The 
first  row  is  associated  with  feedback  from  output  one  which 
is  the  bank  angle  and  the  second  row  is  associated  with 
sideslip  feedback.  Now,  if  it  can  be  assumed  that  the 
control  wheel  contributes  most  of  the  control  surface 
movements  needed  to  control  bank  angle  and  that  the  rudder 
contributes  what  is  needed  to  control  sideslip,  then  it  may 
be  possible  to  zero  the  principle  diagonal  elements  of  the  K 
matrices  without  affecting  performance.  This  should  be  true 
since  this  just  liminates  bank  angle  feedback  to  the  rudder 
and  sideslip  feedback  to  the  control  wheel.  If  this  is  done 
Equations  (5-116)  and  (5-117)  become: 


0.0 

1.337 

Ko  " 

1.766 

0.0 

o 

• 

o 

i - 

0.3145 

K1  3 

0.4154 

0.0 

(5-118) 

(5-119) 


By  no  means  will  this  eliminate  the  cross  coupling  that 
occurs  between  the  control  surfaces  and  the  outputs  in  tne 
state  space  model.  What  Equations  (5-118)  and  (5-119) 
really  say  is  that  the  system  is  not  allowed  to  use  bank 
ang le  feedback  in  determining  how  much  rudder  is  needed  or 
to  use  sideslip  feedback  in  determining  how  much  control 
wheel  deflection  is  needed.  The  affect  of  using  minimum  K 
matrices  of  this  type  is  given  later  since  this  is  one  of 


114 


the  considerations  used  in  determining  the  robust 
controllers  of  the  next  section. 

Robustness  of  The  Control  Laws 

Although  control  laws  are  designed  at  each  flight 
condition,  it  is  desirable  that  a  single  controller  should 
result  in  satisfactory  performance  at  all  flight  conditions. 
A  robust  controller,  if  feasible,  would  reduce  the  need  for 
gain  scheduling.  Several  schemes  are  tried  in  order  to 
derive  a  robust  controller.  Gain  elements  of  the  Kg  and  Kj 
matrices  for  each  maneuver  are  added  together  and  averaged 
in  an  attempt  to  arrive  at  a  robust  controller  for  that 
maneuver.  This,  however,  results  in  very  poor  performance 
of  the  controller  when  used  in  flight  condition  #3,  the  low 
altitude  low  airspeed  case. 

The  method  that  proves  to  be  the  most  successful  is  to 
use  the  controller  designed  for  flight  condition  #3  and 
apply  it  to  the  other  two  flight  conditions.  In  simulation 
of  these  robust  controllers  the  Kg  and  Kj  matrices  are 
minimized  using  the  techniques  discussed  earlier,  and  a 
computational  time  delay  is  included. 

Robust  Controller ,  Coordinated  Turn .  The  controller 
used  in  performing  the  coordinated  turn  with  computational 
time  delay  included  is  defined  by  Equations  (5-95)  thru 
(5-101).  Using  the  guidelines  discussed  earlier,  the  Kg  and 
Kj  matrices  of  Equations  (5-100)  and  (5-101)  become: 


v;^,--.y, 


1 

r 

1.337 

II 

o 

■x: 

0.0 

1.766 

0 . 0 

1 

j 

i 

r 

*i 

0.0 

0.3145 

K1  ■ 

0.4154 

o 

• 

o 

(5-120) 


(5-121) 


Controller  matrices  of  this  form  are  referred  to  as 
minimum  K.  Figures  5-12a  thru  5-12f  show  the  -esults  of 
this  minimized  controller  for  flight  condition  #3.  These 
figures  are  titled  "Robust  Controller  Cond  #3  Coordinated 
Turn"  since  this  is  the  controller  used  for  a  robustness 
check.  Table  5-12  summarizes  the  results  obtained  from 
Figures  5-12a  thru  5-12f. 

Comparison  of  Table  5-12  with  Table  5-6  and  Figures 
5-12a  thru  5-12f  with  Figures  5-6a  thru  5-6f  shows  the 
effect  of  minimizing  the  K  matrices.  From  this  it  is  seen 
that  little  has  changed  and  definitely  nothing  is  lost  in 
reducing  the  K  matrices. 

The  next  step  ir.  determining  robustness  of  this  control 
.law  is  to  apply  it  to  the  other  two  flight  conditions.  This 
is  done  and  che  results  for  flight  condition  #2  are  shown  in 
Figures  5-13a  thru  5 — 1 .3 f  and  summarized  in  Table  5-13. 
Results  for  flight  condition  #1  are  shown  in  Figures  5-14a 
thru  5-14f  and  Table  5-14. 
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FIGURE  5-12c :  Robust  Controller  Flight  Condition  #3 

Coordinated  Turn  (Output  2:  Sideslip) 
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FIGURE  5-12d:  Robust  Controller  Flight  Condition  #3 

Coordinated  Turn  (Control  1:  Rudder) 
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FIGURE  5-13c: 
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Robust  Controller  Flight  Condition  #2 
Coordinated  Turn  (Output  2:  Sideslip) 
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FIGURE  5— 13d : 


Robust  Controller  Flight  Condition  #2 
Coordinated  Turn  (Control  1:  Rudder) 
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FIGURE  5-14c: 
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Table  5-12 

Simulation  Results  For 

Robust  Controller  Condition  #3  Coordinated  Turn 


P 


Table  5-13 

Simulation  Results  For 

Robust  Controller  Flight  Condition  #2  Coordinated  Turn 


I nput/ 
Output 

Peak 

Value 

Final 

Value 

tp(sec) 

ts (sec) 

<t>  (deg) 

30.217 

30.0 

7.4 

5.4 

0  (deg) 

-0.008 

-0.005 

1.0 

15.6 

r  (deg/sec) 

1.38 

1.19 

4.2 

6.0 

4r (deg) 

-■0,46 

— 

— 

— 

4w(<3eg) 

26.0 

— 

— 

— 

Overall,  the  robust  controller  is  able  to  perform  the 
coordinated  turn  with  relatively  good  results.  The  only 
thing  that  might  at  first  appear  unacceptable  is  the  rudder 
displacement  for  flight  condition  #2  shown  in  Figure  5-13d. 
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Input/ 

Output 

Peak 

Value 

Final 

Value 

tp(sec) 

ts (sec) 

■»  t 

<1 

r( 

</>  (deg) 

30.563 

30.0 

7.0 

5.2 

k'< 

7 

/«, 

0  (deg) 

i 

o 

• 

o 

00 

o 

• 

O 

CD 

10.0 

10.6 

‘  '<1 

•'1 

r  (deg/sec) 

4.20 

3.93 

5.5 

8.0 

M 

&r (deg) 

-5.0 

— 

— 

— 

4w(deg) 

30.0 

— 

— 

— 

-3 
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ft 
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At  first  glance  there  appear  to  be  large  osc i 1 lat ions  in  the 
rudder  while  the  maneuver  is  being  performed.  Closer 
examination  shows  these  oscillations  to  be  on  the  order  of 
0.1  degree.  This  small  amount  of  rudder  oscillation  can 
hardly  be  considered  unacceptable. 


Table  5-14 

Simulation  Results  For 

Robust  Controller  Flight  Condition  #1  Coordinated  Turn 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec) 

ts (sec) 

<t>  (deg) 

30.228 

30.0 

7.2 

5.2 

0  (deg) 

-0.0104 

-0.006 

6.8 

14.8 

r  (deg/sec) 

1.35 

1.24 

4.6 

6.0 

fcr  (deg) 

-0.6 

— 

— 

— 

£w(deg) 

2  9.0 

— 

— 

Robust  Controller ,  S idesllp  Using  the  same  approach 
as  for  the  coordinated  turn,  a  robust  controller  is  tested 
for  the  sideslip  maneuver.  Equations  (5-102)  thru  (5-108) 
define  the  controller  used  for  flight  condition  #3  with  time 
delay  when  performing  the  sideslip  maneuver.  The  K  matrices 
are  again  minimized  and  the  simulation  performed.  This 
time,  however,  unacceptable  oscillations  of  the  control 
wheel  occur r .  Thus  in  order  to  keep  the  effect  of  changing 
the  K  matrices  to  a  minimum  Equations,  (5-107)  and  (5-108) 
are  changed  to 
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y. 


II 

o 

0.0 

1.049 

2.486 

-0.9278 

0.0 

0.3497 

K1  = 

0.8285 

-0.3093 

(5-122) 


(5-123) 


Therefore  the  robust  controller  for  the  sideslip 
maneuver  is  defined  by  Equations  (5-102)  thru  (5-104)  along 
•  :.th  Equations  (5-122)  and  (5-123). 

Applying  this  robust  controller  to  flight  condition  #3 
produces  the  results  shown  in  Figures  5-15a  thru  5-15d  and 
Table  5-15.  Results  for  flight  condition  #2  are  given  in 
Figures  5-16a  thru  5-16d  along  with  Table  5-16  while  results 
for  flight  condition  #1  a  e  shown  is  Figures  5-17a  thru 
5-17d  and  Table  5-17. 


Table  5-15 

Simulation  Results  For 

Robust  Controller  Flight  Condition  #3  Sideslip 


Input/ 

Output 

Peak 

Value 

Final 

Value 

tp(sec) 

tg (sec) 

4>  (deg) 

-0.4774 

0.0 

8.4 

— 

0  (deg) 

5.088 

5.0 

16.0 

8.8 

Sr  (deg) 

6.0 

— 

— 

— 

5w(deg) 

28.0 

— 

— 

— 

Again  comparison  of  Table  5-15  with  5-8  along  with 
Figures  5-15a  thru  5-15d  with  Figures  5-8a  thru  5-8d  shows 
that  there  is  little  effect  when  minimizing  the  K  matrices. 
Overall,  the  robust  controller  is  able  to  perform  the 
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sideslip  maneuver  with  fair  results.  However  some  increases 
in  settling  time  are  noted. 


Table  5-16 

Simulation  Results  For 

Robust  Controller  Flight  Condition  #2  Sideslip 


III  1  "  L 

Input/ 

Output 

Peak 

Value 

F  inal 
Value 

t  (sec) 

tg (sec) 

<t>  (deg) 

-0  ,7199 

0.0 

8 . 6 

— 

0  (deg) 

5.06 

5.0 

16.0 

11.6 

3r  (deg) 

8.5 

— 

— 

— 

«w(c<09) 

50.0 

—  — 

— 

— 

Table  5-17 

Simulation  Results  For 
Robust  Controller  Flight  Condition  #1  Side 


Input/ 

Output 

Peak 

Value 

Final 

Value 

tp (sec) 

ts (sec) 

<t>  (deg) 

-0.37993 

0.0 

8.4 

— 

0  (deg) 

5.005 

5.0 

16.0 

8.8 

5r (deq) 

9.0 

— 

— 

— 

5w(deg) 

47.0 

— 

— 

— 
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FIGURE  5-17a :  Robust  Controller  Flight  Condition  #1 

Sideslip  (Command  Input  Vector) 


FIGURE  5-17b: 


I  !  ME .  SECONDS 


Robust  Controller  Flight  Condition  #1, 
Sideslip  (Output  1:  Roll  Angle;  Output 
2:  Sideslip) 
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i obus t  Controlle r ,  Normal  Climb  The  controller  used 
to  perform  the  normal  climb  for  flight  condition  #3  with 
computational  time  delay  is  given  by  Equations  (5-25) 
thru(5-31).  The  robust  controller  used  to  perform  the  climb 
maneuver  for  a’*  fight  conditions  comes  from  this  control 
law  except  .ad  Kj  of  Equations  (5-3Q)  and  (5-31)  are 
changed  to 

-0.9558  0.0 

Kn  =  (5-126) 

0.0  2.957 

-0.3823  0.0 

K,  *  (5-125) 

0.0  1.183 

•  - 

This  represents  the  minimum  form  for  the  K  matrices  of 
this  maneuver.  Notice  that  the  first  element  in  the  second 
row  of  both  Kq  and  K,  have  been  set  to  zero.  This 
eliminates  the  velocity  feedback  from  affecting  the  elevator 
input.  Thus  the  elevator  only  receives  input  from  the 
feedback  of  the  flight  path  angle. 

Figures  5-18a  thru  5-18h  and  Table  5-18  show  the 
results  of  this  robust  controller  when  used  to  perform  the 
normal  climb  for  flight  condition  #3.  Figures  5-19a  thru 
5-19h  and  Table  5-19  show  the  results  for  flight  condition 
#2  while  Figures  5-20a  thru  5-20h  and  Table  5-20  snow  the 
results  obtained  for  flight  condition  #1. 

Comparisons  of  Figures  5-18a  thru  5-18h  with  Figures 
5-9a  thru  5-9h  and  Table  5-18  with  Table  5-9  shows  the 
effect  of  reducing  the  Kq  and  Kj  matrices  to  Equations 
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(5-70)  and  (5-71),  respectively.  Overall,  the  robust 
controller  perforins  well  over  the  three  flight  conditions  if 
the  settling  times  of  flight  condition  #2  and  flight 
condition  #1  are  not  considered  to  be  excessive. 


Table  5-18 

Simulation  Results  For 

Robust  Controller  Flight  Condition  #3  Normal  Climb 


Input/ 

Output 


7  (deg) 
u  (ft/sec) 
Se (deg) 
5t(%RPM) 


Peak 

F  inal 

t  ( sec 

Value 

Value 

P 

5.433 

4.86 

8.1 

0.31 

o 

. 

o 

4.8 

ts (sec) 


12.9 


-3.5 

24.0 


Table  5-19 

Simulation  Results  For 

Robust  Controller  Flight  Condition  #2  Normal  Climb 


ts (sec) 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec 

7  (deg) 

1.759 

1.49 

8.4 

u  (fc/sec) 

-0.375 

0.0 

4.2 

5e  (deg) 

-1.0 

— 

— 

5-j.  ( %R?M) 

30 

— 

— 

15.3 


'7  ‘  '  k  '  > 

-  *  -  ►  b  * 


COMflflW  •  NPUT  VtCTOR 


er  Flight  Condition 
tate  1:  Altitude) 


(  INCHES  ) 


FIGURE  5-18f:  Robust  Controller  Flight  Condition  #3, 

Normal  Climb  (State  7:  1st  Body  Bending 
Mode) 
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Table  5-20 

Simulation  Results  For 

Robust  Controller  Flight  Condition  #1  Normal  Climb 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec) 

tg (sec) 

7  (deg) 

1.786 

1.54 

8.4 

15.0 

u  (ft/sec) 

-0.247 

— 

3 . 6 

— 

Se  (deg) 

-1.25 

— 

— 

— 

5t (%RPM) 

20.5 

— 

— 

— 

Robust  Controller ,  Pitch  Pointing  Pitch  pointing 
proved  to  be  a  difficult  maneuver  to  perform  from  a 
robustness  standpoint.  In  fact  a  robust  controller  that 
produces  satisfactory  results  for  all  three  flight 
conditions  cannot  be  found.  However,  another  purpose  of 

i 

this  section  is  to  show  the  effect  of  minimizing  the  KQ  and 
Kj  matrices  used  to  perform  maneuvers  for  flight  condition 
#3.  Using  the  controller  given  by  Equations  (5-109)  thru 
(5-115)  with  KQ  and  Kj  changed  to  Equations  (5-126)  and 
(5-127)  which  represent  the  minimum  form,  the  maneuver  is 
again  simulated.  Figures  5-21a. thru  5-21h  and  Table  5-21 
show  the  results  of  this  minimum  K  controller.  These 
results  again  show  that  selective  elements  of  the  Kg  and  K 
matrices  can  be  set  to  zero  without  affecting  simulation 
r  asults . 
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RLTITUDE  FT.  i  VELOCITY  f  T /SEC  £  COfIMRND  INPUT  VECTOR 

-2. CO  -t.i>0  -l.CO  -0 . 1)0  O.CO  O.bO  O  O.CO  0.60  1.00  2.40 


TIME.  SECONDS 


JRE  5-21a:  Flight  Condition  #3,  Pitch  Pointing  with  Time 
Delay,  Minimum  K  {Command  Input  Vector) 


O.CO  2.CC 


FIGURE  5-21b: 


6.  CO  8.00  ! 0 , 00 

TINE.  SECONDS 


Flight  Condition  #3,  Pitch  Pointing  With  Time 
Delay,,  Minimum  K  (Output  1:  Altitude;  Output 
2 :  Veloc ity ) 
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ELEVATOR  (DEC)  t  THRUST  (  XRPfl ) 


Flight  Condition  #3,  Pitch  Pointing  With  Time 
Delay,  Minimum  K  (Control  1:  Elevator; 

Control  2:  Thrust) 


FIGURE  5-2  Id  : 


FIGURE  5-21es 


MflE.  SECONDS 


Flight  Condition  #3,  Pitch  Pointing  With  Time 
Delay,  Minimum  K  (Control  3:  Speed  Brakes) 


^  v-v- 

so:- 


FIGURE  5—21 f : 


Flight  Condition  #3,  Pitch  Pointing  With  Time 
Delay,  Minimum  K  (State  7:  1st  Body  Bending 
Mode) 
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Table  5-21 

Simulation  Results  For  Pitch  Pointing 
With  Time  Delay  Flight  Condition  #3  Minimum  K 


Input/ 

Output 

Peak 

Value 

F  inal 
Value 

tp ( sec) 

ts (sec) 

h  (ft) 

0.598 

0.0 

4  ,  6 

— 

9  (deg) 

4.0 

4.0 

16.0 

5.2 

u  (ft/sec) 

-0.899 

0.0 

3.0 

— 

5e  (deg) 

5.0 

— 

— 

— 

5s^deg) 

52  .,0 

— 

— 

— 

<$t(%RPM) 

32.0 

—  — 

“  — 

“  — 

Conclusion 

This  cnapter  presents  the  design  results  for  the 
KC-135A  aircraft.  Design  procedures  along  with  tracker 
control  laws  are  presented  for  each  of  the  three  flight 
conditions  studied.  Simulation  results  in  the  form  of 
Calcomp  plots  for  the  maneuvers  of  flight  condition  #3  are 
also  presented.  Plots  for  the  maneuvers  of  flight 
conditions  #1  and  #2  are  given  in  Appendix  C.  Also  included 
in  the  chapter  are  the  effect  of  computational  time  delay  on 
the  maneuvers  performed  for  flight  condition  #3.  Finally 
robustness  is  checked,  using  what  is  called  minimum  K 


VI.  CONCLUSION  AND  RECOMMENDATIONS 


Thesis  Summary 

This  thesis  is  designed  to  demonstrate  the  use  of 
Professor  Brian  Porter's  multivariable  control  law 
techniques  in  synthesizing  digital  controllers  for  the 
KC-135A. 

Chapter  I  presents  an  introduction  to  the  study  area. 
Chapter  II  gives  a  brief  overview  of  the  control  law  theory 
for  unknown,  known/regular,  known/ irregular ,  and  B*  design 
techniques.  Chapter  III  describes  the  KC-135A  aircraft. 
Chapter  IV  introduces  the  concepts  of  computational  time 
delay  and  presents  the  modifications  made  to  MULTI  to 
incorporate  a  time  delay.  Chapter  V  describes  the  actual 
development  and  analysis  of  the  control  laws  developed  for 
the  KC-135A.  This  includes  presentation  and  discussion  of 
the  effects  of  computational  time  delay  on  control  law 
performance.  Also  included  in  Chapter  V  is  the  development 
of  a  robust  controller  for  each  flight  maneuver,  which  is 
capable  of  performing  a  single  maneuver  over  the  entire 
range  of  flight  conditions.  This  robust  controller  is  of 
the  minimum  K  form. 

Appendix  A  presents  the  development  on  the  aircraft 
equations  of  motion  along  with  tabular  listings  of  stability 
axes  and  body  axes  derivatives  for  the  KC-135A  describing 
each  of  the  three  flight  conditions  studied.  Appendix  B 
gives  an  example  of  the  output  of  the  computer  program  CAT 


used  to  transform  stability  axes  derivatives  to  body  axes. 
Simulation  data  for  each  of  the  maneuvers  performed  for 
flight  conditions  one  and  two  are  presented  in  Appendix  C. 
Appendix  D  describes  a  fix  incorporated  by  this  author  to 
provide  properly  scaled  calcomp  plots  when  using  the 
computer  program  MULTI.  Appendix  E  describes  the  debugging 
of  the  fortran  code  for  MULTI  that  was  necessary  in  order  to 
incorporate  the  computational  time  delay  of  Chapter  IV. 
Appendix  F  presents  a  method  of  transforming  state  equations 
to  the  O-B^  form  of  Chapter  II. 

Conclusions 

This  research  work  demonstrates  that  tracker  control 
laws  can  be  synthesized  for  the  KC-135A  aircraft  using  the 
singular  perturbation  methods  developed  by  Professor  Brian 
Porter.  Tne  controllers  designed  are  able  to  tightly  track 
the  command  inputs  and  at  the  same  time  provide  for  a  high 
degree  of  output  decoupling. 

As  more  constraints  are  placed  on  the  system,  such  as 
control  surface  actuators  and  computational  time  delays,  it 
becomes  increasingly  difficult  to  provide  for  desired  output 
tracking  without  inducing  unacceptable  control  surface 
deflections.  In  almost  all  cases  some  compromise  has  to  be 
made  between  system  performance  and  control  surface 
deflections.  Chapter  V  points  out  the  occurrence  of  these 
compromises  for  the  maneuvers  of  flight  condition  #  3  as 
they  occur. 


Chapter  V  also  introduces  a  reduced  form  for  the 
control  laws  developed.  The  reduced  models  are  called 
minimum  K  designs.  The  results  of  these  minimum  K 
controllers  demonstrate  the  feasibility  of  reducing  the 
required  number  of  feedback  gains. 

Results  of  the  robustness  test  again  show  the 
sensitivity  of  the  control  surfaces  to  parameter  changes. 
This  control  surface  sensitivity,  that  keeps  occurring, 
suggest  the  importance  of  actuator  models.  Simple  first 
order  actuators  are  used  in  this  study  but  higher  order 
models  including  some  sort  of  prefilter  may  be  required. 

The  use  of  MULTI  in  the  development  of  multivariable 
tracker  controllers  produces  relatively  good  design  with  a 
very  moderate  expenditure  of  manpower.  Once  the  user 
becomes  familiar  with  the  use  of  MULTI  and  the  effect  of 
parameter  change  on  design  results,  a  design  can  easily  be 
performed  in  an  afternoon.  The  only  disadvantage  to  this 
method  proved  to  be  the  sensitivity  of  the  control  surface 
deflections,  apparently  caused  by  the  high-qain  nature  of 
the  design. 

Several  days  before  the  due  date  for  this  thesis  an 
error  was  discovered  in  the  fortran  code  for  the  simulation 
portion  of  MULTI.  This  error  has  since  been  corrected.  The 
error  causes  some  of  the  correct  simulation  results  to 
differ  from  those  given  in  Chapter  V.  This  discovery  was 
made  to  late  for  the  results  to  be  included  in  this  thesis. 
However,  some  of  the  changes  are  worth  noting.  First  the 
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settling  times  for  .e  given  maneuvers  increase  slightly, 
approximately  0.2  seconds.  This  however  could  probably  be 


reduced,  if  desired,  by  redesigning  the  controllers. 

Second,  and  most  important,  the  sideslip  manuever  in  Chapter 
V  is  performed  with  a  controller  that  differs  from  the  one 
used  to  perform  the  coordinated  turn.  The  reason  for  this 
change  is  given  in  Chapter  v.  However,  with  the  code 
changed  to  correct  for  the  error,  this  problem  no  longer 
exists  and  the  controller  used  to  perform  the  coordinated 
turn  will  now  work  equally  as  well  for  the  sideslip 
maneuver.  Other  changes  that  occur  are  very  small  and  would 
not  significantly  modify  the  results  given  in  Chapter  V. 

Recommendat ions 

The  recommendations  that  are  presented  in  the  following 
parr  raph  are  a  direct  result  of  problems  encountered  in  the 
completion  of  this  thesis.  These  problems  resulted  in  an 
mdepth  study  of  the  computer  code  for  MULTI,  especially  in 
the  code  involving  the  simulation  and  calcomp  routines. 

Thus,  some  cf  the  recommendations  pertain  to  MULTI  itself. 

Appendix  A  gives  the  state  space  representation  for  the 
lateral  model.  In  the  designs  presented  in  Chapter  V  the 
outputs  are  chosen  as  angles,  specifically  the  roll  angle 
and  the  sideslip.  This  results  in  the  matrix  product  CB 
having  rank  deficiency.  Thus  the  r.eed  for  an  irregular 
design  requiring  a  measurement  matrix.  If  instead  of 
argles,  rates  are  used,  specifically  roll  rate  and  yaw  rate 
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the  matrix  product  CB  would  have  full  rank  and  the 
measurement  matrix  would  not  be  required.  This  design 
however  has  transmission  zeros  in  the  right  half  of  the 
s-plane.  A  design  was  attempted  on  this  model,  however, 
limited  results  were  obtained  since  the  aircraft  proved  to 
be  unstable.  By  increasing  the  sampling  time  to  0.1  seconds 
the  aircraft  became  stable  provided  no  actuators  were  used. 
Thus  the  first  recommendation  is  for  more  research  to  be 
conducted  on  plants  having  unstable  transmission  zeros  in 
order  to  determine  the  feasibility  of  the  design  methods  on 
this  type  of  system. 

An  accurate  means  of  locating  transmission  zeros  for 
plants  not  in  the  O-Bg  form  of  Equations  (2-22)  and  (2-23) 
is  needed.  However,  if  the  plant  can  be  expressed  in  the 
form  of  these  equations,  then  the  transmission  zeros  can  be 
found  by  using  Equation  (2-41)  for  regular  plants  or 
Equation  (2-63)  for  irregular  plants.  Thus  an  option  in 
MULTI  that  would  transform  state  equations  to  the  0  form 
of  Equations  (2-22)  and  (2-23)  would  be  helpful.  This  would 
also  allow  for  easy  implementation  of  Equations  (2-41)  and 
(2-63)  to  determine  the  exact  location  of  the  transmission 
zeros.  Appendix  F  presents  a  simple  method  for  determining 
a  transformation  matrix  that  will  perform  this 
transformation  with  minimum  affect  on  the  state  vector. 

Another  area  that  needs  additional  attention  is  the 
implementation  of  actuator  dynamics.  Current  choices 
include  embedding  them  in  the  A  matrix  as  is  done  in 
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Reference  16,  or  adding  them  to  MULTI  via  Option  #  4. 

MULTI’s  simulation.  Option  #  26,  handles  these  two  methods 
very  differently.  This  stems  from  the  fact  that  MULTI  makes 
a  call  to  ODE  once  during  each  execution  of  the  inner  loop 
of  the  simulation.  Thus,  a  full  analysis  of  the  simulation 
portion  of  MULTI  pertaining  to  these  differences  should  be 
accomplished.  This  indeed  is  necessary  since  the 
conclusions  of  this  thesis  show  that  actuator  dynamics  Dlay 
an  important  part  in  the  design  process. 

When  using  MULTI,  if  an  Option  is  entered  by  mistake, 
it  must  be  allowed  to  finish  since  there  is  no  way  to  exit 
it  without  terminating  the  program.  Reference  17  describes 
another  interactive  program,  TOTAL,  that  allows  for  an  exit 
of  an  unwanted  option  by  inserting  a  $  as  a  response  to  any 
requested  prompt.  This  produces  an  extremely  user  friendly 
atmosphere  that  would  surely  be  beneficial  to  future  MULTI 
users. 

Finally  MULTI  currently  is  restricted  to  100k  of  core 
memory  in  order  to  run  interactively  on  the  Cyber.  Addition 
of  the  recommended  options,  discussed  earlier,  could  cause 
MULTI  to  exceed  this  limit  since  it  is  currently  approaching 
it.  Thus  a  .::-.jor  modification  of  the  program  to  provide  for 
more  efficient  use  of  memory  allocations  is  essential. 
Concentration  of  this  effort  should  be  placed  on  the 
s'nulation  and  Calcomp  routines. 
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APPENDIX  A 


b 


AIRCRAFT  EQUATIONS 
OF  MOTION 


V.*.', 


Introduction 

This  appendix  presents  the  models  of  the  KC-135A 
aircraft  for  the  three  flight  conditions  given  in  Chapter 
III.  The  models  are  developed  using  linearized  equations  of 
motion  and  are  represented  as  continuous  state  space  models. 


Assumptions 

The  following  assumptions  apply  to  the  equations  of 
motion  as  they  appear  in  this  Appendix. 

1.  The  equations  are  assumed  to  have  no  coupling 
between  motion  in  the  lateral  directional  plane  and 
the  longitudinal  plane. 

2.  The  X,  Y,  and  Z  body  axes  lie  in  the  plane  of 
symmetry  and  the  origin  of  the  axis  system  is 
located  at  the  center  of  gravity  of  the  aircraft. 

3.  Aircraft  mass  is  assumed  to  remain  constant  for 
each  flight  condition. 

4.  The  aircraft  is  a  rigid  body. 

5.  The  earth  is  an  inertial  reference. 

6.  The  equations  are  perturbation  equations  and 
represent  the  aircraft  for  small  perturbations 
about  the  trimmed  flight  condition. 

7.  The  flow  is  quasistealy. 

8.  Thrust  acts  only  in  the  X-axis  direction. 
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Equations  of  Motion  (Body  Axes ) 
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Using  the  assumptions  stated  above  the  linearized 
longitudinal  equations  of  motion  are: 
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Where  60  ,  aglj»  and  6^.  represent  displacements  in  elevator, 
speed  brakes,  and  thrust  respectively.  Equation  (a-5)  is 
derived  from  Figure  A-l  using  small  angle  approximation. 


Figure  A-l:  Geometry  of  a  Climb 
Using  the  same  assumptions  the  linearized  lateral 
equations  of  motion  are: 
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Where  (  is  used  to  distinguish  body  axes  from  stability 
axes.  The  equations  used  to  convert  lateral  derivatives  to 
bodv  axes  are: 
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(A-23) 


m 

SS 

s> 

ft 

v'>. 

'\0 

Vs 


-  (Cj  +  Cn  )sinancosaf 


(C^  )b  =  C1  cos  aQ  -  (Cn  -  )  sinagcosa 


0 


.  2 

+  C  sin  a. 

np  0 
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(A-24) 


(A-25) 


y.i’oX'Tv' 


b 


( A  -  2  6  ) 


(C 


CT  cos  cx  -  C  sincr 
Ls  0  ns  0 


(C  )h 


C  cosa.  +  c.  sin  „ 

n0  °  ifl  ° 


!Cnp)fa 


Cnpcos  °0  '  (Cnr  -Clp)sin«0coa«0 

"  cirsin2»o 


(C  )  K  =  C_  cos  a  +  (C,  +  C  )sinancOsan 

nf  b  n  ^  0  lj.  n  p  ^  ^ 

+  C,  nin2a 

P  u 


(Cn  >b  =  Cnrcosa0  +  C1  sina0 
do  o 


(Cy  >b 
Y0 


=  c. 


0 


(CY  >fc 

*  P 


C  cos^  -  C  sinaQ 

3  p  Jr 


<Cy  >b 

J  r 


C  cos  aQ  +  C  sina0 

r  Jp 


!Vb 


=  c. 


(A-27) 

(A-28) 

(A-2  9 ) 
(A-30) 
(A-31 ) 
(A-32 ) 

(A-3  3 ) 
(A-34) 


D imens ionalized  Body  Axes  Equat ions 

Once  the  conversion  to  body  axes  is  made  it  is 
desirable  to  dimensionalize  the  derivatives.  The  equations 
for  dimensionalization  of  the  longitudinal  control 
derivatives  as  given  in  Reference  17  are  as  follows: 


z*  - 

z  Cz 

Zo 

(A-3  5) 

zq 

[ (Z  c)/(2U0)3Cz^ 

(A-3  6 ) 

Zu  = 

(2/u0)cZu 

(A-3  7 ) 
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^  ^  "v  MV 


V  >  % %  r. ' m-  **  -"V  TV  W-' 


m 

k-V 


I 


•5: 

■.V, 

.V, 


&  ' 
V-- 
& 

» 

fa 


Ld  “  L<ci^b 


Lp  =  [  (L  b)/(2U0)]  <cip)b 
L r  =  t  (L  b)/(2U0)  i  <C1  )  b 
L.  ■  L  (C,  ) 


p 


Vb 


Y'Cy  >b 

a 


yp  =  i  < v  bi/(2u0n  (cy  )b 


yr  -  1  (V  b)/(2U0)]  (Cy  )b 


Y«  *  Y(Cy1,b 

0 


(A-52) 
(A-53 ) 
(A-54.) 
(A-55) 
(A-56) 
(A-57) 
(A-58) 
(A-59) 


where  N«  (qSb)  /  (I2Z  )  ,  L=  (qSb)  /  (Ixx)  ,  and  Y=(qSb)/m  and 
q*dynamic  pressure,  s=surface  area  of  wing,  m=mass  of  the 
aircraft,  b=wind  span,  Izz=body  axes  moment  of  inertia  about 
z-axis,  and  Ixx*body  axes  moment  of  inertia  about  x-axis. 


State  Space  Representation 

It  is  desired  to  represent  Equations  (A-l)  thru  (A-9) 
in  state  space  form.  In  doing  so  there  are  a  total  of  nine 
aircraft  states,  two  bending  mode  states,  five  outputs  and 
five  inputs.  The  complete  state  modle  is  obtained  by  using 
Equations  (A-l)  thru  (A-9)  and  dimensional  body  axis 
derivatives  obtained  from  Equations  (A-10)  thru  (A-59)  along 
with  the  following  assumptions: 

1.  X*  =  ZA  =  0 
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>  s.ji  ’j  V-  f  Jf  *-■ 


2.  UQ  =  constant  and  VQ  =  W  =  0 

3  *  *0  ~  P0  =  ^0  “  R0  =  0 

4.  w  =  Uq  and  w  =  Uq<* 

5.  v  =  Ug/3  and  v  =  UQ/3 

The  state  space  model  in  matrix  form  is 


x  *  Ax  +  Bu 


y  =  Cx 


where 


=  [0)8prhguaqee] 


5  [  M  p  r  h  J  j  qee] 

=  [  4  b  6  b  4  ] T 
r  w  e  sb  T 

=  [<0  0  h  e  u]T 
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0 

0 

Y5 

w 

L«, 

Li 

w 

% 

N6 

"w 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

C  is  maneuver  dependent  and 
Appropriate  dimensional  body 
the  prime  notation  using: 


0 

u 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

X  , 

x  . 

X, 

5e 

{8b 

6t 

\ 

K 

5sb 

0 

M’ 

5e 

Ms 

{sb 

0 

0 

0 

0 

0 

0 

discussed  in  Chapter  V. 
es  derivatives  are  changed  to 


x» 

38 

~g  cos8q 

(A-62 ) 

Xq 

= 

Xq  *  U0a0 

(A-63) 

z»' 

= 

{-g/U  0 )  sine0 

(A-64) 

K 

as 

Vuo 

(A-65) 

K 

3S 

<VV  +  1 

(A-66) 

K 

= 

z«/u0 

(A-67) 

= 

Vuo 

(A-6  8 ) 

"# 

= 

Miz, 

(A-69) 

M  * 

u 

+  Mi(Zu/U0> 

(A-70) 

M* 

= 

Ma  +  Mi(Zi/U0) 

(A-71) 

M* 

= 

Mq  +  Mi<<  z,/u0)  +  n 

(A-72) 

'•V’f \"T «lc>V,,\irj< Jr*J*~'K2r+Jr'£r‘*Jfm*J>  vTO  *\ “^V ^N.T'n ^  N.T S.'^S-T'* ^ V *}* ■  *?'-*  *-"\ "* * 


"•  s-*T'  rw  i,',  l 


p 


►:*t 

’  «4 


M  1 

0 

- 

M  +  M  •  ( Z 

0  a 

8/U0> 

(A-73) 

Y* 

= 

g/n0 

(A-7  4 ) 

Y  1 

P 

= 

Vuo 

(A-75) 

yp 

= 

(Vp/u0)  + 

ao 

(A-7  6 ) 

= 

<W  - 

1 

(A-7  7 ) 

= 

Vuo 

(A-7  8 ) 

r  » 

Li  +  ^KZ 

/Xxx>Ni 

/3\»7Q\ 

Li  . 

1 

1  -  <4 

)/(IxxIZ2> 

VA  /  7  ) 

vr  f 

Ni  + 

/Tzz)Li 

u_on\ 

JN  j 

1  -  (1^ 
x  '  xz 

>',<IXXI22> 

O  U  ) 

ud 

= 

(U0)/57.3 

(A-81 ) 

Note  that  i  in  Equations  (A-79)  and  (A-80)  represents 
0,  Pf  r,  Sf  t  and  5W  and  that  i  in  Equations  (A-68),  (A-73) 
and  (A*’78 )  implies  6& ,  <Jr  ,  or  5 w  as  appropriate.  Also 

body  axes  inertias  are  used. 

Actual  transformation  from  non-dimensional  stability 
axes  derivatives  to  dimensional  body  axes  is  accomplished 
using  a  computer  program  created  by  A.  Finley  Barfield 
(Ref.  16)  which  utilizes  the  equations  summarized  in  this 
appendix.  A  sample  run  of  this  program  fcr  flight  condition 
#2  is  shown  in  Appendix  B.  Data  received  from  this  program 
has  the  units  of  radians,  radians  per  second,  and  feet  per 
second.  This  data  is  converted,  by  the  author,  into  units 
of  degrees,  degrees  per  second,  and  feet  per  second.  Body 
axes  derivatives  for  all  three  flight  conditions  studied  are 
listed  in  the  tables  that  follow. 
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Percent  thrust  is  modeled  in  terms  of  thrust  available 
and  is  given  by 

2 

Thrust  Available  ft/sec 

x  =  -  - -  (a-82  ) 

(Mass) 100  %  RPM 

where  m=g/w  is  the  mass  of  the  aircraft  and  thrust  available 
is  measured  in  pounds. 

KC-135A  Aircraft  Models 

The  three  flight  conditions  as  discussed  in  Chapter  III 

are : 

Condition  #1.  High  altitude,  high  speed  cruise  at 

Mach  0.77  and  45,000  feet 

Condition  #2.  Medium  altitude,  heavy  weight  cruise  at 

Mach  0.77  and  28,500  feet 

Condition  #3.  Landing  configuration  at  Mach  0.21  and 

sea  level 

I  artinent  aircraft  data  for  each  of  the  three  flight 
conditions  as  found  in  Reference  10  and  11  is  listed  in 
Table  A-l.  The  nondimensional  stability  axes  derivatives  as 
found  in  Reference  13  are  listed  in  Tables  A-2 ,  A-3 ,  and  A-4 
for  conditions  #1,  #2,  and  #3  respectively.  The  dimensional 
body  axes  derivatives  for  the  A  and  B  matrices  of  Equation 
(A-60)  are  listed  in  Table  A-5,  A-S,  and  A-7  for  conditions 
#1,  #2,  and  #3  respectively.  All  derivatives  not  listed  in 
these  table  are  assumed  o  be  zero. 


Table  A-l 

KC-135A  Aircraft  Data 


Cond it  ion 

#1 

#  2 

#2 

Units 

Altitude 

45,000 

28,500 

sea  level 

ft 

Mach  #  (true) 

0.77 

0.77 

0.21 

— 

Weight 

U) 

««• 

O 

O 

o 

284,000 

130,000 

lbs 

Center  of  Gravity 

32.1 

24.2 

32.1 

%MAC 

q 

124.8 

279.7 

65.9 

lbs/ft2 

S  (wing  area) 

2433 

2433 

2433 

ft2 

b  (wing  span) 

130.83 

130.83 

130.83 

ft 

c  (wing  MAC) 

20.16 

20.16 

20.36 

ft 

UQ  (true) 

745 

745 

745 

ft/sec 

eQ  (body) 

2.4 

2.4 

-0.1 

deg 

aQ  (wing) 

4.4 

4.4 

4.4 

deg 

aQ  (body) 

2.4 

2.4 

2.4 

deg 

* 

*xx 

2,050,000 

2,930,000 

2,050,000 

slug  ft 

xyy 

2,460,000 

4,660,000 

2,460,000 

slug  ft 

Xz2 

4,360,000 

7,480,000 

4,360,000 

slug  ft 

Zxz 

— 

-- 

— 

- — 

Table  A-2 

Non-dimensional  Stability  Axes  Derivatives 
For  Flight  Condition  #1. 


CL 

0.426 

CM 

0.0 

CD 

0.024 

CL 

LU 

0.0 

CM  * 
a 

-6.79 

o.o 

CL 
^ a 

5.329 

CM<, 

-1.1747 

°Da 

0.2417 

CL 

5.1545 

% 

-15.65 

Lq 

°L5 

®e 

0.2114 

CM 

\ 

-0 , 6647 

\ 

0.0 

CL 

-0.3189 

Cm* 

®sb 

0.07259 

CD 

0.0497 

*Sb 

*sb 

% 

-0.223 

C 

n3 

0.166 

% 

-0.762 

Cip 

-0.435 

% 

-0.005 

% 

-0.233 

% 

0.155 

c", 

-0.194 

C^r 

0.428 

% 

0.0315 

% 

r 

-0.113 

\ 

0.264 

Cl*w 

0.0189 

C 

O  A 

6W 

0.00149 

Y5W 

-0.0074 

1 


units  are  radians 
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Table  A-3 

Non-dimensional  Stability  Axes  Derivatives 
For  Flight  Condition  #2. 


CL 

0.426 

CM 

0.0 

CD 

0.024 

ci 

Lu 

0.0 

C»i 

-6.57 

CDU 

0.0 

CLa 

4. 727 

C«. 

-0.8595 

CDa 

0.2143 

S 

4.825 

Cm 

Mq 

-14.65 

S. 

0.1862 

CM 

% 

-0.5988 

CDj 

5e 

0.0 

% 

*sb 

-0.2751 

CM 

M*sb 

0.07639 

CD. 

*sb 

0.04779 

C1 

-0.198 

% 

0.166 

% 

-0.7.2 

C1 

XP 

-0.345 

S 

-0.005 

% 

-0.211 

C1 

r 

0.155 

C_ 

nr 

-0.194 

c*, 

0.428 

c\ 

0.0315 

\ 

-0.113 

\ 

0.264 

Cl*w 

0.0153 

C„ 

n6 

w 

0.00149 

r 

5W 

-0.0074 

units  are  radians 


Table  A-4 

Non-dimensional  Stability  Axes  Derivatives 


For  Flight 

Condition 

#3  . 

0.8108 

CM 

0.0 

CD 

0.0905 

o 

• 

o 

-5.52 

°Du 

0.0 

4.475 

c«„ 

-1.0027 

CD„ 

0.3863 

4.6275 

Cw 

Mq 

-14.05 

0.2222 

-0 . 7105 

C\ 

0.0 

•0.3857 

0.0879 

0.075 

•0.229 

c 

0.132 

Cy/9 

-0.768 

■0.385 

% 

-0.055 

-0.202 

0.248 

% 

-0.186 

0.38 

0.0287 

% 

Jr 

-0.098 

\ 

0.226 

0.0372 

% 

®w 

0.0024 

\ 

-  ; .  0143 

units  are  radians 
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Table  A-5 

Dimensional  A  and  B  Matrix  Coefficients 
For  Flight  Condition  #1. 


X 

u 

-0.0029646 

1/sec 

H 

-4 . 4499 

1/sec2 

0.53477 

ft/sec2 

deg 

LP 

-0.75011 

1/sec 

xi 

-0.53477 

ft/sec 

deg 

Lr 

0.24613 

1/sec 

x; 

-0.56146 

ft/sec2 

deg 

X 

0 . 61895 

i  /  2 

1/sec 

xi 

0.011617 

ft/sec2 

deg 

X 

0.36464 

i  ,  2 

1/sec 

-0.082712 

ft/sec2 

deg 

Nd 

1.42597 

.  .  2 
l/sec 

\ 

0.0495 

ft/sec 
%  RPM 

np 

-0 , 012277 

1/sec 

Mu 

0.010529 

deg 

ft (sec) 

n; 

-0.15052 

l/sec 

-2.7963 

1/sec2 

X 

-0.90838 

1,/sec2 

% 

-0.7537 

1/sec 

N6 

w 

0.20771 

.  /  2 
l/sec 

0.00041339 

1/sec 

y ' 
p 

-0.076917 

l/sec 

X 

-1.64897 

1/sec 

YP 

0.039665 

— 

0.173339 

1/sec 

y; 

-0.99629 

— 

zu 

-0.004883 

deg/ft 

% 

0 . 026647 

l/sec 

K 

-0.54223 

1/sec 

Y5 

w 

-0.00074615 

l/sec 

zi 

0.992966 

— 

Yi 

0 . 04322 

l/sec 

Z  ' 

B 

-0.0018099 

1/sec 

zi 

°sb 

0 . 031962 

l/sec 

X 

-0.021319 

1/sec 

u 

u 

13 . 001745 

ft/sec 

deg 

V 


*yyy*y 


'  r/w 


Table  A-6 

Dimensional  A  and  B  Matrix  Coefficients 
Fight  Condition  #2. 


X 

u 

-0  .  C029479 

1/sec 

L/3 

-6.22193 

.  /  2 
1/sec 

Xa 

0.551527 

ft/sec2 
^  eg  _ 

lp 

-0.904936 

1/sec 

xn 

-0.560066 

f t/sec 
deg 

Lr 

0.382638 

1/sec 

x; 

-0.56146 

f t/sccZ 
deg 

Li, 

0.97056 

1/sec2 

\ 

0.010498 

ft/sec2 

deg 

% 

0.46253 

.  ,  2 
1/sec 

\b 

-0.079813 

ft/sec2 

deg 

N0 

1.87524 

i  /  2 

1/sec 

X{T 

0.034 

ft/sec 
%  RPM 

% 

-0.011695 

1/sec 

M' 

llu 

0.0091035 

deg 

ft (sec) 

n; 

-0.189843 

1/sec 

M' 

a 

-2.40739 

l./sec2 

n; 

r 

-1.18668 

1/sec2 

Mq 

-0.814938 

1/sec 

n: 

sw 

0.02541 

.  ,  2 
1/sec 

0.000441868 

x/sec 

y  ’ 

y/s 

-0.07625 

1/sec 

m: 

6e 

-0.018616 

1/sec2 

YP 

0.0399458 

— 

K 

6sb 

0.0273034 

1/sec2 

y; 

-0.996444 

— 

zu 

-0.0048649 

deg/ft 

ya 

0.02642 

1/sec 

z; 

-0.477268 

1/sec 

y j 

-0  -  0007397 

1/sec 

zi 

0 . 993693 

— 

Yi 

0.0417639 

1/sec 

z'e 

-0.00174889 

1/sec 

K 

6sb 

0.218002 

1/sec 

\ 

-1.758 

1/sec 

% 

13.4555 

ft/sec 

deg 
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Table  A-7 

Dimensional  A  and  D  Matrix  Coefficients 
For  Flight  Condition  #3. 


■0 . 0234233 
0.4228202 


1/sec 

ft/sec 

cleg 


■2 .39756 
■1.11861 


X  ' 

-0.1660307 

ft/sec 

L‘ 

0 . 681707 

deg 

1 

-0 . 5619546 

ft/sec 

deg 

0 .335391 

0.0064484 

f  t/sec 
deg 

L*w 

0.379249 

\b 

-0.063125 

ft/sec 

deg 

Nd 

0 . 588344 

XA 

aT 

0.1268 

ft/sec 
%  RPM 

% 

-0.0852618 

Mu 

0.0318452 

deg 

f  t (sec) 

n; 

-0.238751 

K 

-1.07422 

.  ,  2 
1/sec 

% 

-0.0190339 

-1.09229 

1/sec 

X 

-0.46518 

m; 

-0  00007434 

.  .  2 

1/ sec 

Y0 

-0.129778 

-0.921844 

1/sec2 

y  : 

0.0316453 

K 

4sb 

0.095412 

1/sec2 

Y  ' 
r 

-0.982538 

zu 

-0.0591737 

deg/ft 

0,03819 

K 

-0.778653 

1/sec 

Y‘„ 

-0.0024207 

Zq 

0.966486 

- - 

Yi 

0.1370213 

z'e 

0.000239147 

1/sec 

z ' 

5sb 

0.064588 

\ 

-0.0375145 

1/sec 

u6 

4,101222 

1/sec' 

1/sec 


1/sec4 


1/sec 
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Dody  Bending  (Ref.  10 ) 

One  of  the  assumptions  made  in  the  development  of  the 
equations  of  motion  is  that  the  aircraft  is  a  rigid  body. 

As  a  first  step  in  eliminating  the  rigid  body  assumption, 
this  thesis  models  the  1st  and  2nd  body  bending  modes  for 
the  longitudinal  case.  The  purpose  of  this  model  is  to 
allow  observation  of  the  bending  that  takes  place  while 
performing  a  maneuver,  but  not  to  take  into  account  the 
effects  of  the  bending.  Body  bending  can  affect  sensor 
readings,  such  as  accelorometer s ,  gyros,  ect. ,  if  these 
sensors  are  placed  in  a  position  where  bending  occurs. 
Structural  vibrations  and  body  bending  may  also  have  an 
affect  upon  the  performance  of  the  crew  by  increasing 
fatigue.  The  body  bending  listed  in  this  section  identifies 
crew  station  motions  due  to  structural  elasticity  and 
vibration.  The  structural  bending  mode  equation  is  of  the 
form: 


Ie  +  [(2TwnI)  +  (p0VQCe)]e  +  [w2I  +  pQV2Ce]e 
-  P&V2C  q  *  <,0V\a  -  »0V\  #0 


(A-83) 


where 


=  mass  moment  of  inertia,  slugs-in 
=  elastic  mode  displacement,  inches 
=  viscous  damping 
=  frequency,  radians/second 
=  true  airspeed,  inches/sec 


■  .  -  .  ■  .  •  .  ■  .  .  .TS.-  -  -  ."V  ' 
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C  *  derivative  for  structural  modes 


q  =*  pitch  accelerations,  deg rees/sec^ 
a  *  angle  of  attack,  degrees 

*  elevator  displacement,  degrees 


Using  data  found  in  Reference  10  the  first  and  second 
bending  mode  equations  for  each  flight  condition  become: 

Condition  #1: 

1st  Bending  Mode 

e  +  1 . 313e  +  132»886e  =  -118. 73q  -  148.29or  +  28.496  6e 
2nd  Bending  Mode 

e  +  1 . 1084e  +  543. 60e  =  23.51q  -  10.37a  +  5.518  S& 

Condition  #2 

1st  Bending  Mode 

e  +  1 . 5306e  +  250.593e  =  -268.548q  -  299. 365  **+  48.238  Je 
2nd  Bending  Mode 

e  +  1 . 1959e  +  508. 36e  =  147. 52q  -  l/267«  +  15.759  4e 


Condition  #3 

1st  Bending  Mode 

e  +  1 . 892e  +  131. 34e  =  -62.697q  -  78.306a  +  15.0536e 
2nd  Bending  Mode 

e  +  1 . 225e  +  545. 03e  =  12.428q  -  5.478a  +  2.914*e 


In  adding  the  bending  mode  model  to  the  state  space 
representation  of  the  aircraft,  two  new  states  are  created. 
They  are  e  and  e  where  e  would  represent  the  aircraft 


183 


^r'\^Zr'.r~:.*~*.r ■'r-'.v'-.*r<  •  l^%r-'^*T*T^\T*T'C5' .7:  *  -7'..*  .7^.*” •.'r*.T' \*'-3W'\l 


W. 


h, 

m 

fefi 

V‘0 

to: 

£fJS 

■•a 

r^r*- 

ss 

V’.  V 

tov 

1*W 


to 
ms 

r*vv 

IAO 


V. 'v 


£ 


bonding  modeled  tor  the  crew  station. 

Summary 

This  appendix  summarizes  the  aircraft  equations  of 
motion  used  in  this  thesis.  The  assumptions  used  in 
deriving  the  equations  are  stated,  followed  by  the  equations 
themselves.  Data  for  the  KC-135A  is  in  the  form  of 
non-dimensional  stability  axes  derivatives.  These  are 
changed  to  dimensional  body  axes  derivatives  using  a 
computar  program  described  in  Reference  16  .  Tabular 
listings  of  all  derivatives  are  given  for  each  of  the  three 
flight  conditions  studied.  Also  given  are  the  body  bending 
mode  equations  used  to  depict  bending  of  the  aircraft  at  the 
crew  station, 
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CAIMANS-  ATTACH,  CAT1 ,  IP=T820366 
TFN  IS 
CAT  l 

AT  CY=  001  SN=AFIT 
COMMAND-  CAT l 

-#**»»»  AXIS  iRANSFORMATION  PROBRArt 

ENTER  3TA3IUTY  AXIS  COEFFICIENTS  PGR  TRANSFORMATION 
TO  BODY  AXIS.  TRIM  ALPHA  IS  NEEDED  .'OR  CONVERSION. 

MOMENT  COEFFICIENTS  AND  PIB'FQR'T  COEFFICIENTS  NOT 
RE2UESTED  REMAIN  UNCHAN6ED. 

NOTE:  ALL  COEFFICIENTS  ARE  REQUESTED  WHEN  COMPUTING 
DIMENSIONAL  DERIVATIVES. 

TO  TRANFCQRN  ONLY  LONGITUDINAL  DATA  -  TYPE  LONS 
TO  TRANSFORM  ONLY  LATERAL-DIRECTIONAL  DATA  -  TYPE  LAT 
TO  TRANSFORM  BOTH  LON6  AND  LAT-DXfi  BATA  -  TYPE  BOTH 
KEYWORD  =SQTH 

ARE  Cl MENS I ON AL  BODY  AXIS  DERIVATIVES  REQUIRED  ?  (YES/NOIYES 


!<ll*l(liiSMII«IVVI(riIITITriTTaTITVTTITITITrVmm 


Q  (DYNAMIC  PRESSURE  -  LBS/FT2)  =124.3 

S  (WINS  REFERENCE  AREA  -  FT2)  =2433 

C  (WINS  MEAN  AERODYNAMIC  CORD  -  FT)  =20. 16 

B  (NINE  SPAN  -  FT)  =131.83 

VT  (TRIM  VELOCITY  -  FT/SEC)  =745 

THETA  (PITCH  ANGLE  -  DEBS)  =2.4 

H  I  WEIGHT  -  LBS)  =131109 

INERTIAS  MUST  BE  INPUT  IN  BODY  AXIS. 

I IX  (SLUS-FT2)  =2050000 

IYY  (SLU6-FT2)  =2460000 

111  (SLUG-FT2)  =4360109 

III  (SLUE-FT2)  =0 


HHHHMIHHHHtMIHHHmtHHtHHHHHHti 
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AIRCRAFT  PARAMETERS 

0  (DYNAMIC  PRESSURE  ■■  LBS/FT2)  =  124.803 

S  (KINS  REFERENCE  AREA  -  FT2!  =  2433.00 

C  (KING  MEAN  AERODYNAMIC  CORD  -  FT)  =  20.1600 

I  (KING  SPAN  -  CT)  *  130.830 

VT  (TRIM  VELOCITY  -  FT/SEC)  =  745.000 

THETA  *  2.40000 

N  (HEIGHT  -  LBS)  =  131000. 

Ill  (SLUG-FT2)  *  .203000E+07 

IYY  (SLU6-FT2)  =  .246000E+07 

III  (SLU8-FT2)  =  .436000E+07 

IIZ  (SLUG-FT2)  =  0. 

IS  THE  ENTERED  DATA  CORRECT  ?  (YES/NO! YES 


ALPHA  (DEG)  =2.4 


CL  =.4260 
CLA  (1/DES)  =.093 
CLDE  (1/DEG)  =.003689 
CLDF  (l /DEG)  =-.003367 
CLQ  (1/RAD)  =5.1545 
CLAD  (I/RAD!  *0 
CLU  (1/ (FT/SEC))  =0 
CD  =,0240 

CDA  (1/DEG)  =.004218 
CODE  (1/DE6)  =0 
CDDF  (1/DES)  =.0008674 
CDU  (1/ (FT/SEC) !  =0 

CH  =0 


CMA  (i/DEG)  =-.0205 
CMDE  (1/DEB!  *-.0116 
CMDF  (t/DES)  =.001267 
CMC  U/RAD)  =-13.65 


CHAD  ( 1 /BAD)  =-6.79 


tk^‘  min  *  . 


IT- ' 


V-V. 
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CHU  (1/ (FT/SEC)  I  =0 


LONGITUDINAL  STABILITY  AXIS  COEFFICIENTS 


ALPHA  ~ 

2.40000 

CL  * 

.426000 

CH  » 

0. 

CD  * 

.240000E-01 

CLA  ■- 

.930000E-01 

CHA  « 

-.205000E-01 

Lufl  * 

.421800E-02 

CLDE  » 

.368980E-02 

CMDE  » 

-.U6060E-01 

CODE  > 

0. 

CLDF  * 

-.556700E-02 

CHDF  > 

.126700E-02 

CDDF  « 

.867400E-03 

cm  * 

5.15450 

CHQ  » 

-15.6580 

CLAD  > 

S. 

CHAD  * 

-6.79000 

CLU  > 

0. 

CHU  > 

0. 

CDU  = 

0. 

HHNHHHHHfHHHHIHHHHHtiHHHHHHHttHHHHHHKHHH 

IS  THE  ENTERED  DATA  CORRECT  ?  (YES/NO) YES 


LONGITUDINAL  BODY  AXIS  COEFFICIENTS  (1/RAD! 


■irfSL— ElL 

1 


E1! 

p  v* 

r 


% 

%’ 

S’ 

$ 

K' 

$ 


S; 

Mi 


,  r, 

'S* 


CZ  > 

-.426631 

CX  = 

-.613992E-02 

CZA  > 

-5.37114 

CHA  = 

-1.17353 

CXA  > 

.407432 

CZDE  = 

-.211179 

CXDE  * 

.885101E-02 

CZDF  * 

.316605 

CXDF  « 

-.630117E-01 

CZQ  * 

-5.14998 

CXfi  = 

.215848 

CZAD  > 

0. 

CHAD  > 

-6.78404 

CXAD  * 

0. 

CZU  > 

-.628894 

CMU  > 

•491856E-01 

CXU  * 

-.293671E-01 

LONS 

BODY  AXIS  DIMENSIONAL  DERIVATIVES 

Z  > 

-129542. 

H  > 

0. 

X  > 

-1864.32 

ZA  = 

a  j» •»  naa 

HA  * 

-2.fzfi/ 

XA  = 

30.6425 

ZDE  > 

-15.8825 

NDE  > 

-1.65384 

XDE  » 

.665675 

ZDF  = 

23.8115 

HDF  * 

.180639 

XDF  = 

-4.73904 

ZO  = 

-5.24057 

HQ  > 

-.526903 

Xfi  > 

.219645 

ZAD  = 

0. 

HAD  = 

-.228405 

XAD  * 

0. 

ZU  * 

-.634878E-61 

HU  > 

.164283E-03 

XU  = 

-.296465E-02 

LONG  BODY  AXIS  PRIMED 

DIHENSONAL  DERIVATIVES 

II 

<r 

-.542222 

HA’  = 

-2.79632 

XA'  = 

30.6425 

ZDE’  > 

-.213188E-01 

HDE'  * 

-1.64897 

XDE'  > 

.665675 

ZDF’  » 

.319617E-01 

HDF’  * 

.173339 

XDF’  > 

-4.73904 

ZQ'  » 

.992966 

HQ'  > 

-.753701 

XQ’  « 

-30.9868 

ZU’  > 

-.8521R6E-04 

HU’  > 

.183748E-03 

XU'  > 

-.296465E-02 

ZTHETA’  i 

•  -. 180993E-02 

HTHETA’  = 

.413396E-03 

XTHETA' 

»  -32.1710 

IHHIHHHHIHHHHOtilHIHHItmtHIHIHHfHHHttfHHHHiHH 

CNB  (1/DEB)  >.102897 
CNP  (1/RAD)  >-.005 
CNR  (i/RAD)  >-.194 
CNDR  (1/DEG)  *-.001972 
CNDA  (l/'DEGi  >.000026 
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CNDOT  (1/DEG)  *0 


CROC  (1/DEG) 

CLS  (1/DE6)  *-.00389 
CLP  (1/RAD)  *-.4330 
CUR  (1/RAD)  *.155 
CL DR  (1/DEG)  *.1005497 
CLDA  (1/DE8)  *.0003298 
CLDDT  (1/DE6)  *0 
CLDC  (1/DES)  =0 
CYB  (1/DEG)  =-.013298 
CYP  (1/RAD)  *-.233 
CYR  (1/RAD)  =.428 
CYDR  (1/DES)  =.004607 
CYDA  (l/DE?)  *-.000129 
CYDDT  (1/DE8)  *0 
CYDC  H/DE6)  =0 

LAT'DIR  STABILITY  AXIS  CQEFFICIEM7S 


CND  = 

•289700E-02 

CLB  = 

-.389000E-02 

CYB  * 

-. 132980E-01 

CNP  * 

-.500000E-02 

CLP  = 

-.435000 

CYP  = 

-.233000 

CNR  = 

-.194000 

CLR  = 

.155000 

CYR  * 

.428000 

CNDR  * 

-.197200E-02 

CLDR  * 

.549700E-03 

CYDR  = 

.460700E-02 

CNDA  * 

.260M0E-04 

CLDA  = 

.329800E-03 

CYDA  * 

-.129000E-03 

CNDOT  * 

0. 

CLDDT  * 

0. 

CYDDT  = 

0. 

CNDC  » 

9. 

CLDC  * 

0. 

CYDC  » 

0. 

IHHHCHHIHIWIIIHIHhHItHHHHHHHHilHHHIHKHtHHHiHt 

IS  THE  ENTERED  DATA  CORRECT  ?  (YES/NOYES 


LAT-DIR  BODY  AXIS  COEFFICIENTS 


CN8  = 

.136507 

CLB  * 

-.229636 

CYB  *  -.761919 

CNP  = 

-. 153462E-01 

CLP  * 

-.440833 

CYP  *  -.250718 

CNR  * 

-.188147 

CLR  = 

.144654 

CYR  =  .417868 

CNDR  * 

-.996992E-41 

CLDR  = 

.319410E-01 

CYDR  =  .263962 

CNDA  * 

•227967E-02 

CLDA  * 

.188172E-0! 

CYDA  =•  -.739116E-02 

CNDDT  = 

f. 

CLDDT  * 

0. 

CYDDT  =  0. 
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CNDC  =  0.  CLDC  -  0.  CYDC  =  0. 

#***»****«(*»**»*****»»*<  »»»»*»**'**»«* 


LAT-DIR  SGDY  AXIS  DIMENSIONAL  DERIVATIVES 

NB  = 

1.42597 

LB  = 

-4.4499* 

YB  a 

-57.3031 

NP  * 

-.122772E-01 

LP  = 

-.75*111 

VP  a 

-1.65568 

NR  * 

13*520 

LR  = 

.246128 

YR  * 

2.73949 

NDP  * 

-.9*9384 

LDR  » 

.618955 

YDR  - 

19.8323 

NDA  » 

.2I77I6E-I! 

LDA  = 

.364641 

YDA  * 

-.555881 

NDDT  « 

*. 

LDDT  * 

«. 

YDDT  » 

*. 

NDC  * 

(. 

LDC  * 

». 

YDC  » 

*. 

LAT-DIR  BODY  AXIS  PRIMED 

DIMENSIONAL  DERIVATIVES 

NB’  = 

1.42597 

LB’  * 

-4.44999 

YB’  = 

-,769169E-*1 

NP'  = 

-.122772E-J1 

LP’  = 

-.75*111 

Y?»  = 

•396655E-01 

NR'  = 

15*52* 

LR’  * 

.246128 

YR’  = 

-.996296 

NCR’  = 

-.9*8384 

LDR’  =» 

.618955 

YDR’  a 

.266473E-01 

MDA’  a 

.2*r7*6E-*l 

LDA’  = 

.364641 

YDA’  = 

-.746149E-03 

NDDT’  = 

*. 

LDDT’  = 

*. 

YDDT’  = 

0. 

NDC’  = 

*. 

LDC’  = 

*. 

YDC’  » 

*. 

HHHKHHHtHHiHHOHHHHKHHIHHMIHMHHHHHHHHHHH 

IS  ANOTHER  PROGRAM  RUN  DESIRED  ?  (YES/N01N0 

END  CAT 

•243*1  MAXIMUM  EXECUTICK  FL. 

#.312  CP  SECONDS  EXECUTION  TIME. 
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APPENDIX  C 

S  imulat ion  Results  For 
Flight  Conditions  #1  and  #2 


Table  C-l 

Simulation  Re:  llts  For 
Flight  Condition  #2  Coordinated  Turn 


i 


Input/ 

Output 

Peak 

Value 

Final 

Value 

tp(sec) 

tg (sec) 

<t>  (deg) 

30.48 

30.048 

6.8 

5 . 2 

0  (deg) 

0.0018 

-0.0013 

4.4 

14.0 

r  (deg/sec) 

1.30 

1.19 

4.2 

5.8 

*r  (deg) 

-0.37 

— 

— 

— 

iw(deg) 

22.0 

— 

— 

— . 

Table  c-2 

Simulation  of  Results  For 
Flight  Condition  #2  Sideslip 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec) 

ir 

t  (sec) 

<t>  (deg) 

-0.925 

0.0 

8.4 

— 

0  (deg) 

5.05 

5.0 

16.0 

8.8 

<5r  (deg) 

8.3 

— 

— 

-- 

iw  (deg) 

50 . 0 

— 

— 

Vs  w  v  v  v 

Yk.1  *!V.  mV.  »V.  mV  mV  mV 


Cj 

o 


FIGURE  C-la:  Flight.  Condition  #2,  Coordinated  Turn 
(Command  Input  Vector) 


■v. 


ijL 
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( DECREES/SEX )  m  CONTROL  WHEEL  (DEGREES) 


SURE  C-le 


TIME.  SECONDS 


Flight  Condition  #2,  Coordinated  Turn 
(Control  2:  Control  Wheel) 


- 1 - 

C.CO  2.  CO 

4 ’.CO 

6'.  GO  6.  CO  10.  CO  12.  CO  14.  CO  16. 

|  HUE.  SECONDS 

FIGURE  C-lf: 

Flight 

Condition  #2,  Coordinated  Turn 

(State 

3:  Roll  Rate;  State  4:  Yaw  Rate) 
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ROLL  ANGLE 
2. CO  O.CO 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec) 

P 

ta (sec) 

y  (deg) 

1.73 

1.49 

8.4 

14.7 

u  (ft/sec) 

-0.231 

0.0 

3  .  3 

— 

5e  (deg) 

-1.0 

— 

— 

— 

( %RPM) 

29.0 

— 

— 

— 

Table  C-4 

Simulation  Results  For 
Flight  Condition  #2  Pitch  Pointing 


I nput/ 
Output 

Peak 

Value 

Final 

Value 

t  (sec) 

P 

tg  (sec) 

h  (ft) 

0.147 

0.0 

4 . 4 

— 

$  (deg) 

1.01 

1.0 

16.0 

5.0 

u  (ft/sec) 

-0.192 

0.0 

1 . 6 

— 

SQ  (deg) 

1  .  5 

— 

-- 

— 

4sb(deg) 

18.0 

— 

— 

— 

6_( %RPM) 


48.0 


VELOCir 


FIGURE  C-3 b :  Flight  Condition  #2,  Normal  Climb 
(Output  1:  Flight  Path  Angle; 
Output  2:  Velocity) 


V*  '  \  ",  »  '  »  * 
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ELEVATOR  t  SPEED  BRAKES  ( 
.co  3.00  7.00  n.oo 
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FIGJRE  c-4g:  Flight  Condition  #2,  Pitch  Pointing 
(State  9:  2nd  Body  Bending  Mode) 


Table  C-5 

Simulation  Results  For 
Flight  Condition  #1  Coordinated  Turn 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  (sec) 

t  ( sec 

<P  (deg) 

30.36 

30.0 

6.6 

5.2 

/?  (deg) 

-0.0028 

-0.00237 

7.4 

14.0 

r  (deg/sec) 

1.3 

1.24 

4.2 

5.5 

5r  (deg) 

-0.52 

— 

— 

— 

«w(deg) 

27.8 

— 

— 

— 

Table  C-6 

Simulation  Results  For 
Flight  Condition  #1  Sideslip 


Input/ 

Output 

Peak 

Value 

Final 

Value 

t  ( sec) 

t  (sec) 

<P  (deg) 

-0.665 

0.0 

8.4 

— 

0  (deg) 

5.06 

5.0 

1.6 

8.6 

$r  (deg) 

8.0 

— 

— 

— 

iw(deg) 

48.0 

— 

— 
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FIGURE  C-bd 


Flight  Condition  #1,  Sideslip 

(State  3;  Roll  Rate;  State  4:  Yaw  Rate) 
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Table  C-7 

Simulation  Results  For 
Flight  Condition  #1  Normal  Climb 


Table  c-8 

Simulation  Results  For 
Flight  Condition  #1  Pitch  Pointing 


Input/ 

Output 

Peak 

Value 

Final 

Value 

tp (deg ) 

t  (deg) 

h  (ft) 

-0.16 

0.0 

4.2 

-- 

6  (deg) 

1.01 

1.0 

16.0 

5.2 

u  (ft/sec) 

-0.18 

0.0 

1.9 

— 

Se  (deg) 

1.2 

— 

— 

— 

6sb(deg) 

17.0 

— 

— 

— 

iT(%RPM) 

32.5 

— 

— 

— 

COMMAND  INPUT  VECTOR 
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FIGURE  C-7g:  Flight  Condition  #1,  Normal  Climb 

(State  9:  2nd  Body  Bending  Mode) 


FIGURE  C-7h:  Flight  Condition  #1,  Normal  Climb 
(Sum  of  1st  &  2nd  Bending  Modes) 


APPENDIX  D 


CALCOMP  PLOTS 

Introduction 

As  designs  approached  the  final  stage,  calcomp  plots 
showing  flight  control  system  performance  were  generated 
using  OPTION  #34  and  #35.  Comparision  of  terminal  generated 
plots  and  those  recieved  from  the  calcomp  plotter  showed 
that  plots  whose  y-axis  were  properly  scaled  on  the  terminal 
were  some  times  improperly  scaled  by  the  calcomp  plotter. 
This  was  more  evident  when  attempting  to  plot  curves  that 
had  only  small  deviations  from  zero.  In  some  cases  the  y 
axis  scaling  produced  curves  that  were  hardly 
distinguishable  from  the  zero  line.  See  Figure  D-l  as  an 
example . 


FIGURE  D-l:  Improperly  Scaled  Calcomp  Plot 
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Discussions  with  earlier  users  of  MULTI  verified  that  a 
problem  with  the  calcomp  plots  did  indeed  exist.  However, 
no  one  in  the  past  had  made  any  attempt  to  investigate  the 
problem. 

This  appendix  identifies  a  problem  with  the  calcomp 
plot  routines  found  in  MULTI  and  explains  a  fix  that  is  now 
incorporated  to  eliminate  that  problem. 

Analysis  of  the  Calcomp  Routine 

An  analysis  was  made  of  the  Fortran  code  in  the  plot 
options  of  MULTI  to  determine  why  such  plots  as  the  one 
shown  in  Fig.  D-l  were  being  generated.  When  any  of  the 
options  in  the  thirties  is  called  (i.e.  31,32,33,34,35,35), 
the  first  thing  accomplished  in  MULTI  is  the  generation  and 
storage  of  the  data  to  be  plotted.  All  plot  data  is 
generated  from  the  packed  data  that  is  generated  in  the 
simulation,  OPTION  #26.  This  data  consists  of  a  maximum  of 
101  data  points.  When  the  plot  data  is  generated  it  is 
stored  in  a  matrix  called  PLMAT  which  is  103  rows  by  4 
columns.  At  this  point  rows  102  and  103  are  zero.  Each 
column  of  PLMAT  represents  a  collection  of  data  that  is  sent 
to  the  calcomp  plotter  and  is  used  to  generate  one  curve,  a 
maximum  of  4  curves  can  be  printed  on  each  plot  and  these 
are  explained  in  the  MULTI  Users  Manual,  Reference  4. 

Once  data  has  been  stored  in  PLMAT  the  appropriate 
columns  are  searched  to  determine  the  maximum  and  minimum 
values  for  that  data  being  plotted.  If  two  curves  are  to 
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appear  on  the  calcomp  plot,  then  two  columns  of  PLMAT  are 
searched;  if  three  curves,  then  three  columns  and  so  forth. 
The  maximum  and  minimum  values  as  determined  from  this 
process  are  placed  in  the  number  one  and  two  positions  of  a 
new  vector  called  Y-axis  which  is  103  elements  long.  This 
vector  was  formerly  filled  with  99  data  points  from  the 
second  column  of  PLMAT.  At  this  point  Y-axis  consists  of 
101  points,  where  99  of  them  are  just  fill  points  from 
column  two  of  PLMAT  and  points  102  and  103  are  zero.  Next  a 
call  to  SCALE  is  made  with  the  array  Y-axis  as  the  passed 
data.  SCALE  is  a  subroutine  of  the  calcomp  plotter  that 
determines  the  minimum,  maximum,  and  a  "delta"  increment  for 
a  set  of  points.  The  delta  selected  is  limited  to  1,2, 4,5, 
or  8  times  a  power  of  ten.  SCALE  returns  to  MULTI  the 
Y-axis  data  consisting  of  101  data  points  previously 
generated  plus,  in  position  102,  the  minimum  as  determined 
from  the  call  to  SCALE  and  in  position  103  the  delta  also 
determined  by  SCALE.  Multi  then  assigns  the  102  and  103 
values  of  Y-axis  to  the  102  and  103  rows  of  PLMAT,  in  each 
of  the  columns  to  be  plotted.  Once  this  has  been  completed 
each  column  of  PLMAT  represents  all  data  required  by  the 
calcomp  plotter  to  generate  a  single  curve.  This  includes 
101  data  points,  the  minimum  value,  and  delta  as  determined 
by  SCALE.  The  error  that  caused  plots  as  shown  in  Figure 
D-l  was  generated  by  the  call  to  SCALE.  SCALE  uses  Y-axis 
to  determine  a  minimum  and  a  delta  value.  Y-axis  consisted 
of  99  data  points  that  may  or  may  not  be  plotted  since  they 
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came  from  column  two  of  PLMAT.  This  column  is  plotted  only 
if  two  curves  are  generated  on  the  same  graph.  Since  column 
one  of  PLMAT  would  always  be  used  to  generate  a  curve  on  the 
calcomp  plot,  the  code  was  changed  to  use  this  column 
instead  of  column  two  when  filling  the  data  into  Y-axis. 

This  code  change  eliminates  the  bug  that  existed  in  the 
calcomp  plot  routines. 
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Summary 

An  error  in  MULTI ’ s  Fortran  code  formerly  generated 
calcomp  plot  files  chat  gave  improperly  scaled  calcomp  plots 
under  certain  circumstances.  This  error  was  tracked  down 
and  eliminated. 
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APPEXDIX  E 


Multi  Fortran  Code 


Chapter  4  gives  a  summary  of  the  fortran  code  changes 
necessary  to  incorporate  a  computational  time  delay.  Once 
these  changes  were  made  a  sample  run  of  the  program  was 
accomplished  using  the  examples  given  in  Reference  14. 

These  examples  used  a  computational  time  delay  of  one 
sampling  period.  Comparison  of  the  results  from  the  trial 
run  of  MULTI  and  the  results  given  in  Reference  14  showed 
that  a  problem  initially  existed  with  MULTI.  This  was 
evident  since  MULTI  showed  an  unstable  system  for  the  stable 
system  of  Reference  14.  Thus  it  became  necessary  to  perform 
error  checking  on  the  Fortran  Code  for  MULTI.  As  a  means  of 
doing  this  error  checking,  temporary  print  statements  were 
added  to  the  simulation  portion  of  MULTI  to  print  s(kT), 
r(kT)  and  u(kT)  as  given  by  Equations  (4-1),  (4-2)  and 

(4-3).  A  sample  run  of  multi  was  again  accomplished  using 
the  same  example.  Comparison  of  the  results  of  these  extra 
print  statements  with  ha:.d  calculated  results  showed  that 
the  computer  code  was  performing  the  computational  delay  and 
it's  compensation  correctly  as  given  by  Equations  (4-2)  and 
(4-3).  Further  investigation  revealed  that  the  error 
vector,  e(kT),  used  in  generation  of  the  control  law  was  one 
sampling  period  behind  the  rest  of  the  information  being 
used.  Once  this  small  deviation  was  corrected,  results  from 
f  'LTI  compared  favorably  with  those  given  by  Reference  14. 
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APPENDIX  !•’ 


Transformation  Matrix 


Intrcduct ion 


This  Appendix  presents  a  method,  in  theory  only,  for 
transforming  a  system  of  state  equations  into  the  0-Bo 
form  of  Equations  (2-22)  and  (2-23). 


1  OP 


Transformation  Matr ix 

The  problem  under  consideration  centers  around  the  B 
martix  chat  is  not  in  the  form  of  Equation  (2-22)  and  cannot 
be  put  into  this  form  by  simply  rearrangement  of  the  states. 
An  example  of  this  type  of  matrix  is  given  in  Appendix  A  and 
applies  to  the  lateral  equations  used  in  this  thesis.  Eor 
this  case  the  B  matrix  is  of  the  form 


(D— 1) 


where  B  is  an  nxm  matrix  with  n  equal  to  the  number  of 
states  and  m  equal  to  the  number  of  inputs. 

The  desired  result  is  a  transformation  matrix  T  that 
transforms  Equation  (D-l)  into  the  0-B2  form  of  Equation 
(2-22).  The  proposed  form  for  this  T  matrix  is  shown  in 
Equation  (D-2) . 


0  0 


0  1 
0  0 
0  0 


T23  T24 

1  6  0  i 


(D-2) 
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Notice  the  simple  structure  of  the  T  matrix,  ones  on 
the  principle  diagonal  and  only  two  undetermined  elements,  T 
and  T  .  Notice  also  that  the  B  of  Equations  (D-2) 
contains  the  last  two  rows  of  the  original  B  matrix.  This 
will  always  be  the  case  when  T  is  of  the  form  shown  in 
Equation  (D-2) .  Expanding  Equation  (D-2)  results  in  the 
following  two  equations. 


21 

+  T23331 

+  T24B41 

=  0 

(D- 

-3) 

22 

+  T23B3  2 

+  T24B42 

=  0 

(D“ 

-4) 

Equations  (D-3)  and  (D-4)  represent  a  system  of  two 
equations  and  two  unknowns  (T23  and  T^)  for  which  there  is 
a  unique  solution.  Once  T  is  found  the  new  system  is  given 
by  the  state  equations: 


y  =  C'z 


where 

A'  =  TAT 

c '  »  ct" 

and 


z  *  Tx 


(D-5) 

(D-6) 

(D-7 ) 
(D-8) 

(D-9) 


Notice  that  the  new  state  vector  z  has  all  but  one  of 
the  states  of  the  original  system  and  the  remaining  state  is 
a  linear  combination  of  the  other  states. 
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For  higher  order  systems  the  following  general 
statements  can  be  applied.  Forcing  the  principle  diagonal 
elements  of  T  to  be  exactly  one  results  in  a  system  of 
equations  that  will  always  yield  a  unique  solution  for  the 
unknown  elements  of  T.  The  0-B2  matrix  that  results  from 
the  transformation  will  always  contain  the  bottom  m  rows  of 
the  original  B  matrix.  The  new  state  vector  will  have  n-ui-f 
state  that  are  a  combination  of  the  original  states.  Where 
f  is  equal  to  the  number  of  zero  rows  of  the  original  B 
matrix.  The  remaining  states  will  always  be  the  same  as  the 
original  system. 
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Multivariable  design  technique  developed  by  Professor  Brian 
Porter  of  the  University  of  Salford,  England,  are  used  to  develop 
digital  control  laws  for  the  KC-135A.  Control  laws  are  developed 
for  each  of  three  diverse  flight  conditions.  MULTI,  a  computer 
program  used  in  the  design,  is  modified  to  account  for  computational 
time  delay.  The  effects  of  a  computational  time  delay  on  the 
controllers  developed  are  presented. 

The  controllers  developed,  using  these  techniques,  utilize 
output  feedback  with  proportional  plus  integral  control.  Because  of 
the  structure  of  the  system,  measurement  variables,  in  addition  to 
the  outputs,  are  necessary.  A  reduced  controller  is  modeled  by 
setting  some  of  the  required  feedback  gains  to  zero.  A  comparison 
of  the  results  of  the  reduced  controller  to  the  complete  controller 
is  presented 

A  robustness  controller  is  tested  by  performing  specific 
maneuvers  at  more  than  one  flight  condition.  The  robust  controller 
is  of  the  reduced  form  and  is  required  to  perform  under  the 
constraints  of  the  computational  time  delay. 

This  thesis  concludes  that  actuator  dynamics  play  a  significant 
role  in  the  development  of  control  laws,  and  as  a  result  2nd  order 
or  higher  actuator  models  should  be  used  in  future  studies. 


